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Abstract: For low density gases the validity of the Boltzmann transport equation 
is well established. The central object is the one-particle distribution function, /, 
which in the Boltzmann-Grad limit satisfies the Boltzmann equation. Grad and, 
much refined, Cercignani argue for the existence of this limit on the basis of the 
BBGKY hierarchy for hard spheres. At least for a short kinetic time span, the argu- 
ment can be made mathematically precise following the seminal work of Lanford. In 
this article a corresponding program is undertaken for weakly nonlinear, both dis- 
crete and continuum, wave equations. Our working example is the harmonic lattice 
with a weakly nonquadratic on-site potential. We argue that the role of the Boltz- 
mann /-function is taken over by the Wigner function, which is a very convenient 
device to filter the slow degrees of freedom. The Wigner function, so to speak, labels 
locally the covariances of dynamically almost stationary measures. One route to the 
phonon Boltzmann equation is a Gaussian decoupling, which is based on the fact 
that the purely harmonic dynamics has very good mixing properties. As a further 
approach the expansion in terms of Feynman diagrams is outlined. Both methods 
are extended to the quantized version of the weakly nonlinear wave equation. 

The resulting phonon Boltzmann equation has been hardly studied on a rigor- 
ous level. As one novel contribution we establish that the spatially homogeneous 
stationary solutions are precisely the thermal Wigner functions. For three phonon 
processes such a result requires extra conditions on the dispersion law. We also 
outline the reasoning leading to Fourier's law for heat conduction. 
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1 Goals and Introduction 



Dielectric crystals, as Si and GaAs, have their electronic bands completely filled and 
separated by an energy gap from the conduction band. Therefore electronic energy 
transport is suppressed and the dominant contribution to heat transport is due to 
the vibrations of the atoms around their mechanical equilibrium position. Below 
room temperature these deviations are small, typically only a few percent of the 
lattice constant, hence by necessity weakly anharmonic. As envisioned by R. Peierls 
in 1929 P, the obvious theoretical option is to regard the anharmonicities as a, in 
a certain sense, small perturbation to the perfectly harmonic crystal, which at the 
very end leads to a kinetic description of an interacting "gas of phonons" in terms of 
a nonlinear Boltzmann transport equation. The actual computation of the thermal 
conductivity of dielectric crystals is then based on the phonon Boltzmann equation. 
Through the work of many, for example see [21 El IH E], it has become apparent 
that such a program can be made to work resulting in a reliable prediction over 
a considerable temperature range. Only recently the kinetic description has been 
augmented by molecular dynamics, which numerically solves the classical equations 
of motion, see for example j> . To determine the thermal conductivity one computes 
either the Green-Kubo formula in an equilibrium system at a fixed temperature or 
the average energy flux in the steady state with a temperature difference imposed 
at the boundaries. 

In this note I focus on the step from the weakly anharmonic lattice dynamics 
to the kinetic equation. As an aside, I discuss a few basic properties of the phonon 
Boltzmann equation, mostly to provide some indication on the physics which persists 
on the kinetic level but also to advertise an evolution equation which apparently has 
received little attention. 

If the goal is to compute the thermal conductivity of real crystals, the derivation 
of the Boltzmann equation is considered as a minor issue, where the emphasis varies 
from author to author. Much more relevant is to have reliable information on the 
lattice structure, on the phonon dispersion law, and on the lowest order anharmonic 
elastic constants. Furthermore, on the kinetic level the conductivity is determined 
through the inverse of the linearized collision operator, which cannot be computed 
by hand. Hence suitable approximation schemes had to be developed. I will have 
nothing to say on these topics. 

On a qualitative level kinetic theory provides a rather simple picture for the tem- 
perature dependence of the thermal conductivity, k(T). At "high" temperatures a 
semiclassical approximation suffices, which predicts k(T) = with some temper- 
ature independent coefficient 9^. At "low" temperatures the quantization of lattice 
vibrations must be taken into account. The total number of phonons then equals 
J d 3 k(e^ k '^ kBT — which reflects the freezing of the number of energy carriers as 
T — > 0. On the other hand also momentum nonconserving collisions become rare, 
resulting in a phonon mean free path which diverges as T — > 0. This latter effect 
dominates and yields the prediction n{T) = e e ^ T , 6\ > 0, as T — > 0. Experimentally 
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Figure 1: Thermal conductivity of Si (natural abundance) [7] . 



such a behavior is masked by the finite size of the sample and only over a narrow 
temperature range the exponential increase in 1/T can be seen. A crucial point in 
the experiment is to manufacture a crystal which has no dislocations and is free 
of impurities. Even then, isotope disorder provides an additional mechanism for 
diffusive energy transport, which persists in the harmonic approximation. E.g., for 
Si the natural abundance is 28 Si 92.23%, 29 Si 4.76%, and 30 Si 3.01%, which means 
that the deviation from the perfect constant atomic mass crystal can be considered 
as small. 

To provide an example we reproduce in Figure ^ the thermal conductivity for 
chemically pure and dislocation free Si as measured by Glassbrenner and Slack 
jjj- On the right hand side the importance of the various scattering mechanisms 
is displayed. Above 100°K one notes the classical 1/T-behavior. Below 100°K 
the quantization of phonons becomes relevant. Diffuse boundary scattering reflects 
the size of the probe which is 2 cm long times 0.44 cm as average diameter. The 
umklapp scattering refers to momentum nonconserving collisions, see Section EJ The 
experimental findings are well reproduced by the theory which is based on the 
linearized Boltzmann equation, as will be explained in Section 

In the kinetic theory of gases the central object is the Boltzmann distribution 
function Nf(r,v,t), N the total number of particles, which counts the number of 
gas molecules in the volume element d 3 rd 3 v in the one-particle phase space close 
to r,v at time t. Phonons are not such local objects. In fact, upon specifying the 
complete displacement field, including its velocities, it is not so clear how to extract 
from it the positions and momenta of the particle-like objects called phonons. Most 
likely, for a general displacement field no such procedure can be devised. Still in 
the kinetic limit the mechanical picture becomes precise. As has been recognized 
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for some time jHUHj, the link between a wave field and transport equations allowing 
for a mechanical interpretation is provided by the Wigner function. This approach 
will be followed also in these notes, noting already now that the collision between 
phonons, while they conserve energy and momentum, are otherwise unlike collisions 
between mechanical point particles. 

For the purpose of a better understanding of the validity of the kinetic descrip- 
tion, my guiding principle is to discard all details and to devise the arguably simplest 
of all models, which still displays the same physics. I will even go as far as to ignore 
the obvious fact that atoms deviate in three-space from their equilibrium position. 
Hence I will assume that the displacement field is scalar. The virtue, I hope, is to 
make the derivation of the transport equation maximally transparent. 

We propose to ignore quantization in the first round. One reason is the hope that 
for a classical model techniques different from a hierarchy of correlation functions 
and Feynman diagrams might become available. As a further bonus, we establish 
the link to weakly anharmonic, in general multicomponent, wave equations, which 
are applied in the wave dynamics of the upper ocean, in acoustic turbulence, and 
in other areas ^Oj- In this context the phenomenon of interest is a turbulent state 
maintained through external forcing. Again, kinetic theory is the natural theoretical 
tool to explain and predict properties of the steady state. 

Acknowledgements: I am most grateful to Jani Lukkarinen for many instructive 
discussions and a first reading of the notes. I thank Carlo Cercignani for help towards 
the H-theorem and Eric Carlen for discussions on the Brout-Prigogine equation. 
These notes were first presented as lectures at the workshop "Quantum Dynamics 
and Quantum Transport", Warwick, September 6 - 12, 2004. I am grateful to Gero 
Friesecke for this opportunity. 



2 A real crystal simplified 

We consider the simple cubic lattice Z 3 as the lattice of mechanical equilibrium 
positions of the crystal atoms. The deviations from their equilibrium position are 
denoted by 

q x eR, xeZ 3 , (2.1) 
with the canonically conjugate momenta 

p x eR, xeZ 3 . (2.2) 

We will use units in which the mass m of an atom equals one. For small deviations 
from the equilibrium position we may use the harmonic approximation in lowest 
order. The corresponding potential energy then reads 

t4arm(g) = - ~ V^xQy ■ (2-3) 
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The elastic constants a(x) satisfy 

a(x)=a(-x), \a(x)\ < a e~ ailxl (2.4) 
for suitable oto, ol\ > 0, and 

J2 «(*) = (2.5) 

x& 3 

because of the invariance of the interaction between the crystal atoms under the 
translation q x q x + a. Mechanical stability requires 

a(k) > for k ^ (2.6) 

for the Fourier transform a of a. 

The anharmonicity is assumed to reside only in the on-site potential which we 
divide into a harmonic piece and the rest 

C/site = + Kn(fe)) ■ (2.7) 



Physically, the on-site potential is artificial and it would be more natural to assume 
that the atoms are coupled through a weakly anharmonic pair potential. As we will 
argue below, in the kinetic limit only the collision rate turns out to be modified. 
Thus, for the purpose of deriving the kinetic equation, we might as well stick to the 
somewhat simpler on-site potential. 

The Hamiltonian of the anharmonic lattice system is written as the sum 

H — H + V . (2.8) 

H is the harmonic piece given through 

xGZ 3 x,y£Z s 

uo > 0. The lowest order type of anharmonicity reads 

v = V M V M = ( 2 - 10 ) 

xez s 

with A small. The potential energy £/harm + ^site is then not bounded from below, 
which however will not be visible on the kinetic time scale. If preferred, one could 
add to V an the quartic term X'q x with A' = A 2 /18c<Jq. Then H > and the quartic 
term disappears in the kinetic scaling. For reasons of readability we will set A' = 0. 

We work in the physical space dimension. Whether the kinetic approximation is 
valid in one and two dimensions remains debated. On the other hand only for such 
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low dimensional systems extensive numerical results are available, to which we will 
turn in Section [T71 

The equations of motion are 

Px(t) , 

- ^2 a{y - x)q y (t) - ui*q x (t) - \q x (t)\ x G Z 3 . (2.11) 

We will consider only finite energy solutions. In particular, it is assumed that 
\Px\ ~ * 0) I Qx I ~ * sufficiently fast as \x\ — > oo. In fact, later on there will be the 
need to impose random initial data, which again are assumed to be supported on 
finite energy configurations. As to be explained in great detail, in the kinetic limit 
the average energy diverges suitable linked to the nonlinearity A — > 0. 

We will mostly work in Fourier space and set up the notation. Let T 3 = [— |, |] 3 
be the first Brillouin zone of the dual lattice. For / : Z 3 — > K we use the following 
convention for the Fourier transform, 

?{k) = Y,e~*** h '*fx, fcef. (2.12) 

f(k) extends to a 27r-periodic function on R 3 . The inverse Fourier transform is given 
by 

f x = [ dke i2 ^ x f(k) , (2.13) 

where dk is the 3-dimensional Lebesgue measure. This convention has the advantage 
of maximally avoiding prefactors of 2tt. The dispersion relation for the harmonic 
part Hq is easily computed as 

u(k) = (u 2 + a(k)) 1/2 . (2.14) 

By mechanical stability u(k) > ujq- If ujq > 0, then uj is a real analytic function on 
T 3 . If ujq = 0, uj may still be real analytic, one example being a(k) ~ |/c| 4 for small 
k. In Fourier space the equations of motion become 

p(k, t) , 
-uj(k) 2 q(k,t) 

-X dkxdk 2 S(k - h - k 2 )q(ki,t)q(k 2 ,t) (2.15) 

with k G T 3 . Here 5 is the 5-function on the unit torus, to say, 5(k') carries a point 
mass whenever k' G Z 3 . 



d , x 

d . . 
dt Pxit) 



9^, . 
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It will be convenient to concatenate q x and p x into a single complex- valued field. 
We set 

a(k) = -^(y^{k)q{k)+i^=p(k)) (2-16) 

with the inverse 

q(k) = l=-±=( a (k) + a(-kY) , p(k) = - a(k) + a(-k)*) . (2.17) 

The a-field evolves as 

d f 

—a(k,t) = —iu}[k)a(k,t)—i\ / dk\dk 2 b{k — k\ — k%) 
at Jf6 

(8u(k)u(k 1 )u(k 2 ))~ 1/2 (a(k 1 , t) + a(-k u t)*)(a(k 2 , t) + a(-k 2 , t)*) . (2.18) 
In particular for A = 0, 

a(k,t) = e- luj(k)t a(k) . (2.19) 

For real crystals the a-field would be vector-valued for two reasons: the displace- 
ments are in M 3 and the unit cell contains usually more than one atom. Correspond- 
ingly io then becomes a fc-dependent matrix. Furthermore by translation invariance 
the potential energy of the crystal depends only on the differences q y — q x . As long 
as the interest is merely in the derivation of the Boltzmann equation such extra 
features can be ignored. 

If one simplifies anyhow, the reader may wonder why we do not switch to the 
continuum wave equation. In our context a natural option would be the Klein- 
Gordon equation with a weak quadratic nonlinearity, 

d 2 

—(p(x,t) = A(j)(x,t)~ujQ^(x,t)-\(j)(x,t) 2 , xeR 3 . (2.20) 

Another possibility would be the standard wave equation with a cubic nonlinearity 

d 2 

— <f){x,t) = A(f>(x,t) - A0(x,t) 3 , xeR 3 . (2.21) 

We will discuss continuum equations in Section |H1 from which it will become clear 
that the underlying lattice structure plays a crucial role. 

Having agreed upon the basic model (|2.11j) of our enterprise, we have reached a 
point of bifurcation. Physically we should quantize ()2.8j) . together with ()2.9j) . ()2.1()|) . 
according to the standard rules and then investigate the effects small anharmonic- 
ities. On the other hand it seems to be worthwhile not to hurry so much and to 
explore the classical model, which has an interesting structure of its own. In addi- 
tion there could be help from the theory of nonlinear wave equations, which would 
put our claims on firmer ground. Thus in Sections E3 to El we treat the derivation 
of the Boltzmann equation for the classical model. The same program is repeated 
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for the quantized crystal in Sections El and El] with the approach through Feynman 
diagrams explained in Section ^2 In Section [S] we address wave turbulence which is 
concerned with continuum wave equations such as (|2.2(J|) and (|2.21|) . Sections [7| and 
[T2l study properties of Boltzmann equation, in particular the H-theorem. The non- 
linear part concludes with a discussion of the thermal conductivity In the final part 
of our notes we investigate the harmonic crystal with random isotope substitution. 

3 Local stationarity, Wigner function 

The kinetic theory of dilute gases relies on the scale separation between typical 
interatomic distances and the mean free path. As a consequence locally, in regions 
of linear size much larger than atomic distances and much smaller than the mean 
free path, the statistics of particles is Poisson in a good approximation. If f(r,v,t) 
denotes the Boltzmann distribution function at time t, then close to r the particles 
are uniformly distributed with density p(r) = J d 3 vf(r,v,t) and their velocities 
are independent with common distribution f(r,v,t)/p(r). The Poisson distribution 
is singled out from all other conceivable distributions, because it is stationary in 
time with respect to the free gas dynamics, translation invariant in space, and has a 
strictly positive entropy per unit volume. In fact, there are no other such probability 
measures [T2j . 

To transcribe this kinetic picture to weakly interacting phonons, as building 
blocks we need, on the phase space of microscopic configurations {q x ,Px, x G Z 3 }, 
probability measures which are invariant under the free dynamics generated by H , 
stationary under lattice shifts, and have a strictly positive entropy per unit volume. 
The obvious candidates are Gaussian measures with zero mean. By translation 
invariance, their covariance reads 

(q x q y ) = Q(x - y) , (p x p y ) = P(x - y) , (q x p y ) = C(x - y) . (3.1) 

Let Q, P, C denote the corresponding Fourier transforms. Then Q(k) > 0, Q(k) = 
Q(-Jfe), P(k) > 0, P(k) = P(-Jfe), C(k) = C(-k)% and |Cf < QP. Stationarity in 
time yields in addition the relations 

P = u 2 Q, C(x) = -C(-x) , i.e. C(k) = -C(-k) . (3.2) 

Such properties are more concisely expressed through the a-field. Stationarity in 
space-time is equivalent to 

( a (k)) = 0, (a(k)a(k')) = 0, {a{k)*a{k')) = W{k)5{k - k') . (3.3) 

W(k) > and, by convention, W(k) is a 27r-periodic function on R 3 . Inserting the 
definition (J2.16j) and comparing with (|3.2j) results in 

h W (k) + W(-k)) = UuQ(k) + -P(k)) = ooQ(k) = -P(k) , 

A Z, UJ UJ 

^(W{k) - W{-k)) = id{k) . (3.4) 
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The meaning of the covariance W is grasped better by considering expectations 
of some physical quantities. Let us first study the local energy H x , for which we 
equally divide the potential energy between the two elastically coupled sites. Then 

H x = \vl + l^lql + - a ( x ~ v)QxQv (3-5) 

and 

H =Y,H X . (3.6) 

xez 3 

Clearly 

(H x ) = [ dku{k)W{k) . (3.7) 



T 3 



To probe further, we study the flow of energy out of a big box A C Z 3 . Setting 
H A = J2x&A Hx one finds 

J t HK= \Y^ Yl a ( X - y)(-^Py + QyPx) ■ (3.8) 
zeA yez 3 \A 

Since the coupling is not only nearest neighbor, the division into local currents is 
somewhat arbitrary. To be specific, let us choose as one face of A the coordinate 
plane {x, x l = 0}. Then, in the limit A — > oo, the one-component of the energy 
current becomes 

A = o a ( x ~ J/)( - 9(* 1 .o,o)P» + %P(x\0,0)) ■ (3.9) 

x 1 <0y 1 >l 

Upon averaging, using ()3.4|) and (|3.9jl . 

(j e ) = J_ f dkVa{k)W{k) = ^- [ dk{uVu){k)W{k) . (3.10) 



4-7T J T 3 27T Jj3 

Thus it is natural to regard W as number density in wave number space. ujW is the 
energy density and (27r) _1 V 'uj(uW) is the energy current density. Note that if W is 
even, the total energy current vanishes. 

A further important quantity is the entropy per unit volume, which on general 
grounds is defined as the logarithm of the phase space volume at prescribed values 
of the "macrovariables" , see Appendix 118.21 for further discussion. Here we use an 
equivalent short-cut and compute the Gibbs entropy of the Gaussian measure with 
covariance given through W . To do so let us choose the periodic box [1, £} 3 , £ integer, 
and consider the finite volume analogue of the Gaussian measure from ()3.3|) . Then k 
takes the discrete values k e (^ _1 [1, ...,£]) 3 . Let pc be the corresponding probability 
density. As usual, the entropy of pc is given through 

S e = - //VWogp G = Yl (log^(fc) + log7T + l) (3.11) 

R<3 fce^-Mi,-/]) 3 
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and thus the entropy per unit volume by 

lim r 3 S e = [ dk( log W (k) + log n + l) . (3.12) 

&-»00 J T 3 

The next step is to construct, out of the Gaussian measures introduced in (j3.3J) . 
Gaussian measures which have a slow variation in physical space Z 3 and which are 
locally stationary. For this purpose we give ourselves the local power spectrum 
W(r,k) > 0, r G M 3 , which vanishes rapidly as |r| — > oo, and introduce 

Q(r,x) — / dkW(r, k)uo(k)~ l cos(27rfc ■ x) , 

P{r,x) = / cf&VF(r, &)a;(&) cos(27rA; • x) , 

C{r,x)= dkW{r,k)sm{27ck-x), (3.13) 

a; G Z 3 , by which we define the family (-) ' e of Gaussian measures through 

(g,) G ' £ = o, ( Px } G ' £ = o, 

(q x q x ') G ' £ = Q(e(x + x')/2, x-x')+ 0{e) , 
(p x p x ') G ' e = P(e(x + x')/2, x-x') + 0{e) , 
(q x p x >) G ' £ = C(e(x + x')/2, x - x') + 0(e) . (3.14) 

The error of order e has to be allowed so to ensure a positive definite covariance 
matrix. 

This family has two important properties. 

(i) Relative to the reference point r/e, r G M 3 , the measure becomes stationary in 
the limit e —>■ 0. This is the property of local stationarity. 

(ii) For two distinct reference points r and r', r ^ r', the local distributions become 
independent in the limit e — > as can be inferred from 

lim { (qir/ej+xQlr/ei+x'Qlr'/ej+yQlr'/ej+y> )°' £ 

-(q[r/e\+ X q[r/el+ X ') G ' £ (q[r'/el+ y q[r'/s\+y') G ' £ } =0, (3.15) 

with [_-J denoting integer part, since Q(r,x) — > as |a;| — > oo . The analogous 
property holds for the remaining covariances. Thus under the Gaussian measure 
(■) G ' £ two macroscopically far apart regions are statistically independent. 

The construction ()3.14|) is computationally not so flexible and it is more con- 
venient to invert the order. Thus the primary object is a family (-) G,e of Gaussian 
measures (non-Gaussian measures to be discussed further on). They have mean zero 
and a local covariance, which is almost time stationary and slowly varying in space. 
These conditions are most easily imposed through the lattice analogue of the local 
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power spectrum W expressed in terms of a-field, compare with (|3.3J) . Firstly we 
require 

(a(k)) G > £ = , (a(k)a(k')) G > £ = . (3.16) 
The local a*a spectrum is defined through 

W\x, k) = 2~ 3 [ d7]e i27TX - v (a(k - r]/2)*a(k + r//2)) G ' £ . (3.17) 

(a(k — r]/2)*a(k + r]/2)) G,£ is T 3 -periodic in k and (2T) 3 -periodic in rj. Therefore 
W(x, k) as inverse Fourier transform with respect to 77 is T 3 -periodic in k and lives 
on the half-integer lattice (Z/2) 3 with respect to x. 

We rescale the lattice to have lattice spacing e through the substitution x = e~ 1 y, 
y G (eZ/2) 3 , and obtain the rescaled local power spectrum 

W £ (y, k) = {e/2f [ dr]e l27Ty - v (a(k - eri/2)*a(k + er//2)) G ' e . (3.18) 

Then, denoting |_-J £ as modulo e, one requires 

\imW £ ([r\ £ ,k) = W(r,k) (3.19) 

pointwise. If (-) G,£ is defined through (J3.14)) . then W of (|3.19j) agrees with the one 
in f)3. 13j) . W £ (y, k) is normalized as 

V f dkW £ (y,k)= [ dk(a(k)*a(k)) G ' £ . (3.20) 

ye(eZ/2)3 JTi ^3 

The condition that the limit in ()3.19|) exists thus implies that the average phonon 
number increases as e~ 3 , equivalently the average total energy increases as 

/ dku(k)(a(k)*a(k)} G ' £ = (H } G ^ = 0(e- 3 ) . (3.21) 

()3.18|) has a familiar touch. Recall that for a quantum wave function ip on 
physical space M 3 the Wigner function is defined by 

W £ (x,k)=[ dr}t lx ^{k - er]/2)*ii( y k + er]/2) (3.22) 

with x, k G 1R 3 and ip the Fourier transform of ip. e is the semiclassical parameter, 
e — ► in the semiclassical limit. The main difference to (|3.18j) is that for the 
semiclassical limit usually one considers a sequence ip £ of wave functions, while in 
(J3.18J) one has a sequence of probability measures over the wave field and its time 
derivative. Because of this obvious analogy we call (J3.18|) the Wigner function, more 
properly the one-point Wigner function. The n-point Wigner function is understood 
as the n-th moment of a* a. 
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For a family (-) £ of general measures on phase space one defines the one-point 
Wigner function 

W £ (y, k) = {e/2f [ drie i2ny - v (a(k - eri/2)*a(k + er]/2)) £ , (3.23) 

i.e. through ()3.18|) with (-) G ' e replaced by (-) e . The rescaled two-point Wigner 
function becomes 

W £ {yi, h, y 2 , k 2 ) = {e/2f \ drjxd^ exp[i27n/i ■ r) X + i2ny 2 ■ r] 2 ] 

J(2J/e) 6 

(a(h - ern/2)* a(h + er ll /2)a{k 2 - er] 2 /2)*a{k 2 + er] 2 /2)) £ , (3.24) 
and similarly for higher-point Wigner functions. We require f!3.19|) and 

m n 

lim(na(%rn a ^)) £ = ( 3 - 25 ) 

3=1 i=l 

whenever m ^ n. The condition of statistical independence of far apart regions then 
reads 

lim{W e (LriJ ej A;i, L^Je, fc 2 ) ~ W £ ([n\ £ , k x )W £ {\r 2 \ £ , k 2 )} = (3.26) 

for r x r 2 , which in the context of low density gases is known as assumption of 
molecular chaos. Since ()3.26|) is a law of large numbers, it implies that 

n 

lim W £ ([n\ £ , k h ..., [r n \ £ , kn) = T[W([r 3 \ £ , kj) (3.27) 

£^0 ■*■ ■*■ 

i=i 

whenever the family {r 1; . . . ,r n } is free of double points. 

There is no reason that (-) £ becomes locally stationary as e — ■> 0. Still the 
condition of local stationarity can be expressed through the limiting behavior of 
multi-point Wigner functions. For example, in the case of the two-point function 
the condition would read 

lim W £ ( [r\ £ , k u [r\ £ , k 2 ) = W(r, k x )W(r, k 2 ) 

£^0 

+5{h + k 2 ) [ dr ] W(r,k 1 +ri/2)W(r,k 2 + ri/2). (3.28) 

For a sequence (-) G ' £ of Gaussian measures satisfying ()3.19j) the identity ()3.28|) holds 
by construction. 
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4 Kinetic limit 



As initial measures for ()2.11|) we adopt the scale of Gaussian measures (-) G,£ satis- 
fying ()3.16|) - ()3.19|) . The time-evolved measure at time t is denoted by {-) t . Let us 
first consider the harmonic lattice dynamics, A = 0. Then by linearity, (-) t is again 
Gaussian. Since the deviations from stationarity are on the spatial scale e~ x and 
since there is a finite speed of propagation, one has to wait for times of order e~H to 
observe appreciable changes of the Wigner function, which defines the kinetic time 
scale e~H. On that scale one has 

^- t (a(k-er ] /2ya(k + e7 1 /2)) t/£ , 

= -is" 1 (u{k + e?7/2) - u(k - en/2)) (a{k - er]/2)*a(k + eq/2)) t/e . (4.1) 

Taking the limit e — > one obtains 

d — — 

j-W{r], k, t) = -iVuu(k) ■ r]W(r], k, t) (4.2) 

and, upon inverting the Fourier transform, the limit Wigner function is the solution 
of the transport equation 

%-W(r, k, t) = ~Vu(k) ■ V r W(r, k, t) . (4.3) 
at 2tt 

Thus in the kinetic limit, e — > 0, we can think of the phonon counting function W as 
arising from a gas of independent particles, the phonons, with kinetic energy uj{k). 
Detailed proofs for the validity of the free transport equation ()4.3|) are given by 
Mielke [THj. He allows for rather general deterministic initial data and for harmonic 
lattice dynamics with vector displacements and a general unit cell. 

If one adjusts the strength of collisions in such a way as to have an effect of 
the same order as the transport term, then kinetic theory claims that the locally 
stationary state imposed at t — retains its structure in the course of time. Of 
course, the time-evolved measure (-)t/e is no longer exactly Gaussian. But for small 
e and on a local scale it does remain so in a good approximation. As crucial difference 
to (J4.3|) the evolution equation will contain a collision term taking account of the 
anharmonicities. As to be shown in the following section, the cubic term is of the 
right strength if one substitutes 

A ~* v^A (4.4) 

with A fixed and independent of e. Then the stabilizing quartic term has the strength 
A' = (A 2 /18u;g)£:, which is indeed small compared to the cubic term. The Wigner 
function at the kinetic time t is given through 

W £ (y, k, t) = e 3 [ drie i2ny - r > '(a(k - er]/2)*a{k + erj/2)) t/e . (4.5) 



14 



Figure 2: Three phonon collisions. 



It is expected that the limit e — > exists, 

\imW £ ([r\ E ,k,t) = W(r,k,t), (4.6) 

e— >0 

and the limit phonon counting function W is the solution of a Boltzmann-like equa- 
tion. Its derivation will be explained in the section to follow, but let us state the 
result already now, 

JW, k, t) + -L Vu(k) ■ V r W(r, k, t) (4.7) 
at 2tt 

= 7 / dk 1 dk 2 (u(k)u(k 1 )u(k 2 ))- 1 {25(u(k) + - u(k 2 ))5(k + k t - k 2 ) 

(W(r, fc x , t) W(r, fc 2 , t) + W(r, k, t)W(r, k 2 , t) - W(r, k, t)W(r, h, t)) (I) 
+ 6(u(k) - u(ki) - uj(k 2 ))5(k -h- k 2 ) 

(W(r, ki, t)W(r, k 2 , t) - W(r, k, t)W(r, k u t) - W(r, k, t)W(r, k 2 , t)) } (II) 
with 7 the strength of the collision term, 

7 = \ A 2 ■ (4-8) 

5 is the torus 5-function, see the explanation below Equation (|2.15|) . As shorthand 
the collision operator is denoted by C(W). 

Dynamically the terms (I) and (II) can be viewed as given in Figure El In 
(I) the phonon with wave vector k collides with a phonon with wave vector k\ in 
order to merge into a phonon with wave vector k 2 . The loss term is the term 
proportional to W(k), hence 2W{k)(W{k 2 ) — W(ki)), and the gain term is the 
remainder, i.e. 2W (ki)W (k 2 ) . Note that the gain term has a definite sign while the 
loss term takes both signs. Correspondingly in (II) the phonon with wave vector 
k splits into two phonons with wave vector k\ and k 2 . The gain term is again 
W{ki)W{k 2 ) and the loss term is — W(k)(W(k 1 ) + W(k 2 )). The precise way of how 
the phonon distribution functions appear in (J4.7j) does not seem to have a mechanical 
interpretation in terms of colliding point particles. As can be seen from the 8- 
functions in the collision operator, in both collision processes energy is conserved, 
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while momentum is conserved only modulo integers. E.g. for term (I) the 5-function 
yields the constraint 



k + k 1 = k 2 + n, neZ 3 , k,k 1 ,k 2 eT 3 . (4.9) 

In case n = one speaks of a normal process while in case n 7^ of an umklapp 
process. 

The rates appearing in (I) and (II) come out of the computation to be presented 
in Section However, their relative strength 1/2 is required in order for energy to 
be locally conserved. 

Note that the Boltzmann equation preserves the positivity of W. Obviously the 
free streaming term has this property. If W first hits at some point k, W(r, k, t) = 
0, then dW(r,k,t)/dt > 0, d/dt denoting the total time derivative, due to the 
positive gain term and vanishing loss term of the collision operator. Hence at that 
point W cannot turn negative. 

This seems to be a good moment to return to the issue of a potential energy 
which depends only on the differences in the displacements, as would be the case 
for a real crystal. Then u = and V of (|2.1()j) is replaced by 



1 3 

^^^B^-?,) 3 , (4.10) 

x&Z 3 a=l 

ei, &2, e 3 the standard basis of Z 3 , which expressed in terms of the a- fields becomes 

V 3 = A- / dk 1 dk 2 dk 3 5(k 1 + k 2 + k 3 ) 
3 

J] {^(k^iexp^nk*} - l)(a(kj) + a(-%)*)} . (4.11) 
i=i 

Compared to V of (|2.10|) . only the weight in Fourier space has changed. Thus the 
Boltzmann equation remains as in (|4.7jl provided the collision rate (a;(fc)a;(fc 1 )u;(/c2)) _1 
is replaced by 

3 3 

H^uikj^iexp^Trkfl h = k. (4.12) 

j=l a =l 

Close to the origin this collision rate is more singular than the one in (|4.7j) . But 
the general properties of the Boltzmann equation, as to be discussed in Section 
remain in force. 

We hurried a little bit to write down the Boltzmann equation. So the reader 
might wonder why we claim that on the kinetic time scale local stationarity is 
maintained. The point is that the free dynamics, generated by Hq, does not tolerate 
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deviations from local stationarity as long as the free dynamics is given some time act. 
Such a property has been studied in considerable detail by Dobrushin et al. [13] , for 
recent improvements see ^Hj- Roughly speaking, they consider initial measures (-) £ 
on phase space for which (q x ) £ = = (p x ) £ and for which the Wigner function W £ of 
()3.23j) has a limit as in (j3.19|) . In addition they require that under (-) £ spatial regions 
separated by a distance £ with 1 <C £ <C e^ 1 are in essence statistically independent. 
(■) £ is non-Gaussian, in general. This initial state is evolved under the dynamics 
generated by H . Then, for times t where 1 <C t <C e -1 , the Wigner function does 
not change. However locally the oscillators adjust such that the measure becomes to 
a very good approximation Gaussian and satisfies the conditions (j3.16j) and (|3.19|) . 
Thus for times which are short on the kinetic scale the harmonic lattice dynamics 
forces local stationarity. 



5 Conditions on the dispersion relation 

Our discussion seems to indicate that the kinetic description holds independently 
of the particular form of the (short ranged) harmonic interaction potential, in other 
words independently of the (analytic) dispersion relation. As far as the convergence 
to locally stationary Gaussian measures is concerned, this impression is well sup- 
ported |14j . However, for three-phonon collision processes it cannot be taken for 
granted to have a non-vanishing collision operator. If one sets 

E q (k) = uj{k) + uj{q) — uu(k + q) , (5.1) 

clearly conservation of energy can be satisfied only if 

E q (k) = (5.2) 

admits solutions when considered as a function on T 6 . For nearest neighbor coupling 
only, to say a(e) = —1 for |e| = 1, a(0) = 6, and a(x) = otherwise, the dispersion 
relation reads 

^ 1 /2 

co(k) = (w 2 + 2^ (1 -cos(2vrF))) , k = (k\k 2 ,k 3 ). (5.3) 

3=1 

As shown in Appendix 118.11 for this choice E q (k) > ujq/2 > 0. The physically most 
obvious model does not admit three-phonon collisions. 

From this perspective one might wonder whether ()5.2j) can be satisfied at all. An 
example which can be checked still by hand is given by 

3 

w(Jfe) = lu + 2 (! - cos(2ttF)) . (5.4) 

3=1 
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It corresponds to the harmonic couplings 



(6 + w )(g(i,o,o) - %f , -«(?(2,o,o) - 9o) 2 , -(9(i,i,o) - Qo) 2 , 



(5.5) 



all others determined by isotropy and translation invariance. Note that the next 
nearest neighbor couplings are destabilizing. Clearly E (0) = ujq while for q = 
(1/4)(1, 1, 1), k = (1/8)(1, 1, 1) one has E q (k) = lu - 6(^2 - 1) < provided u is 
not too large. 

To have a nonvanishing collision operator we require 



There seems to be no simple sufficient criterion on u, which would ensure ([5.6)1 . 
Numerically one plots E q (k) for random choices for q to find out whether E q (k) 
takes negative values which then implies (|5.6)1 . 

Observe that E q (0) = ujo for all q G T 3 . If ujq > 0, by continuity there is then 
a neighborhood Ao of defined through A = {£; G T 3 , E q (k) > for all q G T 3 } 
and G A . If W(r, k) is supported in A for every r, then C(W) = 0. The free 
flow leaves this set of Ws invariant and therefore such W's evolve merely by free 
streaming. In general, there will other components of T 3 where no collision partner is 
available. In addition, there can be components Ai, A2, . . • , A m such that if k G Aj 
it will remain so under any sequence of collisions. Then in each Aj the system 
equilibrates in the long time limit, but in general the equilibration temperature will 
differ from component to component. For this reason we introduce the notion that 
k G T 3 \ {0} is linked by a collision to q G T 3 \ {0} if E q (k) = 0. Clearly, linkage is 
symmetric. 

Ergodicity Condition (E): For every k, k' G T 3 \ {0} there is a finite sequence of 
collisions such that k is linked to k' . 

In particular for every k 7^ there is at least one collision partner q 7^ such that 
E q {k) = 0. A necessary condition for ergodicity to hold is u = 0. If in ()5.4|) we set 
Uq = 0, then ergodicity is satisfied with one intermediate collision, as can be seen 
from an explicit computation. 

There is a further condition related to the issue of existence of solutions of the 
Boltzmann equation ()4.7|) . If ||W||oo denotes the sup- norm, the collision operator 
can be trivially estimated as 




(5.6) 




(5.7) 



provided ujq > 0. By standard methods of kinetic theory, if 




(5.8) 
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then the Boltzmann equation (J4.7)) has a unique bounded solution for < t < to 
with suitable t . If (|5.8|) does not hold, resp. if u = 0, to establish the existence of 
solutions local in time would require more efforts. 

For the dispersion relation (J5.4|) the condition (J5.8j) is satisfied. In general, for 
(J5.8|) to hold E q (k) has to be a Morse function uniformly in q. To see why, assume 
that E q is not a Morse function, and try to locate points q where the integral in 
(|5.8|) diverges. On the level set {k, E q {k) = 0} one must have 

V k E q {k)=0, (5.9) 

which can be solved locally to yield q = q{k). Thus 

E q{k) (k) = (5.10) 

must have solutions. Secondly the Hessian of E q (k) must have at least one vanishing 
eigenvalue which leads to the condition 

det(Hess E q (k)) = at q = q(k) . (5.11) 

The surfaces in T 3 defined through the level zero sets in ()5.10|) and ()5.11|) will 
generically intersect along a curve. Thus we must be prepared that the integral in 
()5.8|) diverges along a curve in T 3 . Again, no simple sufficient criterion is available 
to ensure (|5.8j) . 



6 Derivation of the phonon Boltzmann equation 
(classical model) 

The textbook derivation of the Boltzmann equation starts from the quantized theory 
as to be discussed in Section El and uses the Fermi golden rule to compute the 
transition rate, see f2] for a particularly lucid discussion. While such a procedure 
yields the correct rates, it provides little theoretical insight why the Fermi golden 
rule would be applicable in such a field theoretical context. Of course, the best of all 
possibilities would be to have a mathematically rigorous derivation. We are far from 
such a goal at present. Instead we offer in this section a derivation based on the 
concept of local stationarity through which higher order correlations can be suitably 
decoupled, see [TH] for a similar argument in the case of a weakly interacting Fermi 
gas on the lattice. Physically, this seems to me the most transparent procedure, 
admittedly with the disadvantage that the approximate local stationarity cannot be 
checked directly. A more systematic approach uses Feynman diagrams, as will be 
explained in Section ITT1 

To properly argue for the validity of the Boltzmann equation (J4.7|) . it is conve- 
nient to work in atomic units for a while. We give ourselves the Wigner function 
W(r, k) > and assume that the initial measure, (-)o, is Gaussian satisfying ()3.16|) 
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and (|3.19J) . The average with respect to the measure at time t is denoted by (-) t . 
We introduce the shorthands 

a(k, 1) = a(k)* , a(k, -1) = a(k) , (6.1) 

and 

(f)(k, fa, k 2 ) = \(8cu(k)u(fa)uj(k 2 ))- 1/2 . (6.2) 
Then the equations of motion (J2.18)) can be written in the more compact form 



d 
dl 



a(k,a) = iauu(k)a(k,a) + iy/ea \_\ I dfadk 2 (f){k,fa,k 2 ) 

x5(— ok + aiki + a 2 k 2 )a(fa, ai)a(k 2 , c 2 ) , o = ±1 . (6.3) 



The two-point function satisfies 

j t (a(p)*a(q)) t = i{u(p) - w(g))(a(p)*a(g)) t + y/£F(q,p,t) (6.4) 



with 



F(q,p,t)—i ) / dfadk 2 ((p(p, ki, k 2 )5(—p + aifa + cr 2 k 2 )(a(fa, ai) 

xa(k 2 ,a 2 )a(q)) t - <j>(q, fa, k 2 )5{q + aifa + cr 2 )(a(p)*a(fa, ai)a{fa, cr 2 ))tj • (6.5) 
We need a second iteration, which we write in integrated form as 

F(q,p,t) = F hom (q,p,t) + y/e dsG(q,p,t - s, s) . (6.6) 

Jo 

The homogeneous term in ()6.6)) reads 

F hora {q,p,t) = i [ dfadfae^^ 2 ^ 

x (4>{p, fa, k 2 )5(-p + aifa + a 2 k 2 )e- itw{q) (a(fa, a x )a{k 2 , a 2 )a(q)) 

-<f)(q, fa, k 2 )5(q + axfa + a 2 )e ituj{p) (a(p)*a(fa, cn)a(k 2 , a 2 )) ) 

= 0, (6.7) 

since in the initial measure odd moments vanish. We conclude that 
d /** 

— (a(p)*a(q)) t = i(uj(p)-uj(q))(a(p)*a(q))t + sj dsG(q,p,t - s, s) . (6.8) 

Following ()3.23j) one switches to Wigner function variables and sets 

W £ ( V , k, t) = £ 3 (a(k - er)/2)*a(k + e V /2)) t/e . (6.9) 
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Then 

-9 

—W £ {r t , k, t) = ie-\uj(k - erj/2) - u(k + erj/2))W £ (r), k, t) 
at 

rt/e 

+e 3 dsG{k + erj/2,k-er]/2,e~ 1 t-s,s). (6.10) 
Jo 

Assuming that W e (r), k, t) converges to W(rj,k,t) as e — > 0, the remaining task is 
to establish that the inhomogeneous term on the right converges to the collision 
operator (|4.7|) acting on W(rj, k,t). 

(1) Local stationarity , Gaussian approximation. The integrand of the inhomogeneous 
term in ()6.6|) is given by 

G(q,p,t,s) = S~] / dk 1 dk 2 V] / dhdl 2 (j)(p, k±, k 2 )4>(q,h,l 2 ) (6-11) 

x (5{-p + (Tikx + o 2 k 2 )5{q + nh + r 2 l 2 )e- ituj{q) 

+5(q + o~\k\ + a 2 k 2 )5(-p + nh + T 2 l 2 )e itw( pA 
xe u ( . M ki)+«Mk2)) {a ( ku ai)a ( k2} a2)a{ i u Tl)a{ i 2: T2))g 

—2 2, / dk\dk 2 2, I dl\dl 2 

x U(p, k h k 2 )6(-p + o-ih + a 2 k 2 )e~ itw{q) (a(q)a(k 2 , a 2 )a(h, n)a(l 2 , t 2 )) s 
-<f){q, ki, k 2 )5{q + o~xk\ + a 2 k 2 )e ituj{p) (a(pYa(k 2 , a 2 )a{h, n)a(l 2 , t 2 )) ( 



As our basic assumption, in the kinetic scaling regime, the average (-) s at the ar- 
guments in question is in a good approximation a locally stationary measure. If so, 
the averages appearing in (|6.11|) can be substituted by Gaussian pairings. Using the 
shorthand k for a(k,a), the approximation amounts to 

(hk-zhl^s = (kih) s {k 2 l 2 ) s + (hh) s (k 2 h) s + (kik 2 ) s {hl 2 ) s , 
(qk 2 hh)s = (qh)s(k 2 l 2 ) s + (qh) s (k 2 h) s + (qk 2 ) s (hl 2 ) s , (6.12) 

and correspondingly for p. By symmetry, upon inserting in ()6.1H) . the first two 
terms on the right are identical and will yield the gain and loss term, respectively. 
The third pairing is subleading and vanishes as e — > 0. Accordingly we set 

G = Gg a i n + Gloss + C su b • (6.13) 

(2) Gain and loss term. In G ga i n we change to Wigner fucntion variables as 

ki = k' — eairj' /2 , h = k' + e<j\r] ' /2 , n = —o~i , 

k 2 = k" - eo- 2 7]"/2 , l 2 = k" + ea 2 r]"/2 , r 2 = -a 2 . (6.14) 
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The e-dependence of <fi can be ignored and the 77-integration is extended to M 3 , since 
by asumption W £ (r], k,t) has a good decay in 77. The phases have to be expanded 
to first order in e. Then 

Gg a i n (£: + er]/2, k - erj/2, t, s) 

= 2e 6 V / dk'dk" [ ^d I ,>(M',fc") 2 e if(fflu(fc '- £ffl,i ' /2)+ * (t "- OT2,, " /2) » 

x (s{-k + £(77/2) + <7i(A/ - + <7 2 (A;" - ea 2 r/'/2)) 

S(k + £(77/2) - cr^A;' + - a 2 (k" + ea 2 r]" y 2 )) e -^( fc W2) 

+5(A; + £(r//2) + ax (A;' - £ C r 1 7 7 72) + <7 2 (fc" - ea 2 r)" /2)) 

S(-k + 5(77/2) - cri(fc' + e<Ti7772) - <r 2 (k" + ea 2 r]" /2))e lt ^ k - £r>/2) ^ 

{a(k' - ea 1 r]'/2, ai)a(k' + ea 1 r]'/2, — <7i)) s 
x(a(k" - ea 2 r]"/2, a 2 )a(k" + ea 2 rj"/2, -a 2 )) s 

= 2£- 3 V / dk'dk" [ drfdrf<]>(k, k', fc'')2 e «(^«(fc')+««(fc")) 



Cl,<72 = 

-ite(Vw(fc)(r?/2)+Vw(fc')(»772)+Va;(fe")(J7"/2)) 



(e-^W5(A; - (nib' - a 2 k") + e^ (fc) 5(A; + tr^' + a 2 k")) 
5(rj -r]' - rj")W £ (r]', k', s)W £ (r]", k", s) 

= 2e- 3 V / cta'dA;" / drfdrf'<j>{k, k', jfe") 2 ( e «(-«(*)+^«(*')+^(*")) + c . c .) 
CTliCT2=±1 iT6 iR6 

,J(jfc _ ^Jfe' _ (T2 jfe") e -*te(V«(fc)(T?/2)+Va;(*')(^/2)+V«(fc")(»?"/2)) 

5(77 - 77' - r]")W £ (r]', k', s)W £ (r]", k", s) , (6.15) 

where in the last step the c.c. term arises through replacing the sum over a±, a 2 by 
the sum over — <7i, — o" 2 - 

In Gioss we change to Wigner function variables as 



dk 3 5(q - k 3 )(k 3 h) s (k 2 l 2 ) s , 

T 3 

h — k' — er]'/2 , k 3 = k' + £t//2 , n = 1 , 

k 2 = k" - ea 2V " /2 , l 2 = k" + ea 2 rj"/2 , r 2 = -a 2 , (6.16) 



and 



dfc 3 5(p - k 3 ){k 3 l 1 ) s {k 2 l 2 ) s , 

T 3 

k 3 = k' - ei]'/2 , h = k' + er)'/2 , n = -1 , 

k 2 = k" - ea 2 r]"/2 , h = k" + ea 2 r]"/2 , r 2 = -<7 2 . (6.17) 
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Then 

G loss (k + erj/2, k - erj/2, t, s) 



= -4 £ 6 V / dkrdk'dk" I dr ] 'drj"(f){k 1 ,k',k")(t)(k,k 1 ,k")a 1 

e it(o-ia)(fex)+<72w(fe // )) e ite(-Vw(fe / )(»//2)-Va)(fc")(^72)) 

(e-iMk') 6 (_ aih + k >_ _ a 2 k" - e{r)"/2)) 

x5(-k + + o x k x + a 2 k" - e(r]"/2))5(k + e(r]/2) — k' — e{rf /2)) 

_^{k') 5 ^_ aiki _ k > _ £ ^y 2 ) _ a2k „ _ 

x5(k + + a x k x + a 2 k" - e(rj" /2))5(k - e(rj/2) — k' + e(r]'/2)) 

W%r)',k',s)W £ (r)",k",s). (6.18) 

We integrate over k' and neglect the shift of order e in the //-argument of W £ . In 
the second summand a±, a 2 is substituted by —a\, — cr 2 with the result 

G\ oss (k + ^/2, k - erj/2, t, s) 

= -Ae- 3 V] [ dhdk" [ dV^ / 0(^^i^ // ) 2 (e^ W ~ CTl£4,(fcl) ~ CT2a,(fe '' )) + c.c.) 

x5(?7 - V - rf')W e {r{, k, s)W £ (t]", k", s) . (6.19) 

By assumption the Wigner function is varying on the kinetic scale. Thus the 
remaining time integration for G ga i n and Gi OS s is of the generic form 

lim J dug{u) dsie^-^ + e- iuj{t ~ s) )f(es,e(e-H - s)) 

lim / dujg{uo)e^ 1 2 I ds cos(us/e)f(t — s, s) 
£ ^°J Jo 

= 2tt / dug(cu)5((j)f{t,0), (6.20) 



where g{oS) is some smooth test function of rapid decay. 

Combining ()6.15|) . ()6.19|) . ()6.20|) and upon noting that the convolution becomes 
multiplication in position space, one concludes that 

lim / ' ds(e/2f I dT]^ lrie ' v (G £ in (k + er]/2, k - £ V /2, (e~H -s),s) 

£ ^°Jo i(2T/ £ )3 v 

+G £ loss (k + er)/2, k - erj/2, (e" 1 * - s), s)) 
= A 2 — / dkidk 2 (uj(k)uj(ki)uj(k 2 )) ~ l 5(u — V\LO\ — o^ 2 ) 

8(k - axh - a 2 k 2 ) (W(r, h, t)W(r, k 2 , t) - 2a 1 W(r, k, t)W(r, k 2 , t)) , (6.21) 
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which agrees with the collision term (|4.7)l . 

(3) Subleading terms. There are two subleading terms from (|6.12|) . denoted here by 
C su b = G su bi + G su b2- For G su bi we change to Wigner function variables as 

k x = k! - ea x rf /2 , k 2 — k' + ea x rf j2 , a 2 = —a x , 

l x = k" - enr]"/2 , l 2 = k" + er x r)" /2 , r 2 = -r x . (6.22) 

Then 

G suhx (k + er] /2, k - erj/2, t, s) 

= V / dk'dk" [ ^W>(^,^,fc , )0(fc^ ,, ^ ,, )e" ite(V " (fc ' ) ' 7 ' +V " (fc){r ' /2)) 

-*M*)«y(_jfe + _ £r /)6(& + - er/') 



+e^ (fe) 5(A; + 5(77/2) - £r/)6(A; - e(ri/2) + £77")) W e (^, fc', s)W £ {rj", k" , s) 

= Ae~ 3 [ dk'dk" [ dri'dri"(f)(k, k', k')<f>(k, k", k") (e" iMfe) + e itu>i - k) ) 
Jt 6 Jr 6 

s ^ e -M^(k'W + ^ { k)( v /2)) 6 ^ _ r] ,_ rj'^w'tf, k 1 , s)W £ (ri", k" , s) . (6.23) 
The remaining time integration is of the generic form 

t/e ft 

dscos{u(0)(£-H- s))f(es) = e _1 / dscos{u(0)s/e)f(t - s) 
Jo 

= cuiO)- 1 ( sin(u;(0)t/e)/(0) + J ds sm{u{0)s/e)f{t - s)) . (6.24) 

The second summand is of order e. The first summand oscillates fastly around zero 
average and thus vanishes by one further integration in time. 

Our argument indicates that o>(0) > is required. If u(0) = 0, then the product 
5{k)u{k)~ 1 is not defined. Whether this is an artifact of the derivation or signals a 
limit in the validity of the kinetic description remains to be understood. 

For G su b2 we change to Wigner function variables as 

dk 3 5(q - k 3 )(k 3 k 2 ) s (l x l 2 ) s , 

k 2 = k'- e V '/2 , k 3 = k' + erf /2 , a 2 = 1 , 

l x = k" - er x rf'l2 , l 2 = k" + 6T 1 r ] "/2 , r 2 = -r x , (6.25) 

and 



dk 3 5(p - k 3 ){k 3 k 2 ) s (l x l 2 ) s , 

T 3 

k 3 = k'- erf 12 , k 2 = k' + erf 12 , a 2 = -1 , 

l x = k" - er x rf'l2 , l 2 = k" + 6T 1 r ] "/2 , r 2 = -t x . (6.26) 
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Then 



G snh2 (k + er]/2, k - er]/2, t, s) 

= -2 V / dhdk'dk" I dri'dr ] "<i)(k,k 1 ,k')<i)(k 1 ,k",k") 

<r 1 ,n=±l JrB 6 
e -ite(V«(fc')(»?'/2)+Vw(fc)(»j/2)) ( j^_ (Tifci _ er li^ ^ e it(aMki)+^{k')-^{k)) 

x5(k + e(rj/2) — k' — e(rj' /2))S(-k + e(rj/2) + a x k x + k' — e(rj'/2)) 

x5(k + e(rj/2) + aifci — k f — e(rf /2))} a^W^rf , k', s)W £ (r]", k", s) 
= -4 V / dhdk'dk" [ dri'dr ] "<i)(k 1 ,k,k)<i)(k 1 ,k",k") 

e ito lW (fci) ^ _ fc / + e ^ 2 ) - £(77'/2)) - 5(fc — k' — £(77/2) + £(7772))) 
(J(<riifei + er]")e- it£V ^ k ' w 5(r] - 77' - ri")aiW e (ri', fc', s)W e (ri", k" , s) . (6.27) 

Integrating over yields the phase u(er]"). If u;(0) > 0, the remaining time integra- 
tion is of order 1. The difference of (^-functions in the large round bracket is of order 
£, when integrated against W £ . Therefore the second sub leading term vanishes as 
e -> 0. 



7 Some properties of the classical phonon Boltz- 
mann equation 

(i) Energy. The energy at position r and time t on the kinetic scale is defined 
through 

e(r,t)= I dku{k)W{r,k,t) . (7.1) 
It satisfies the local conservation law 

|U(r,t)+V-j c M) = 0. (7.2) 
From the transport term one concludes that the energy current is given by 



j e {r,t) = (27T)- 1 / dk(Vu(k))u(k)W(r,k,t). (7.3) 
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The vanishing of the contribution from the collision term can be seen from 
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dk 1 dk 2 dk 3 (uj(k 1 )uj(k 2 )uj(k 3 ))- 1 {25(uj(k 1 ) + uo(k 2 ) - w(k 3 ))6(ki + k 2 - k 3 ) 

00(h) {W(k 2 )W(k 3 ) + W{k x )W{k 3 ) - W(k l )W(k 2 )) + S(u(k{) - u{k 2 ) - u(k 3 )) 
S(kt -k 2 - ^(^(Wik^Wih) - WikJWfa) - W(h)W(k 3 ))} 

= / dkidk 2 dk 3 (u(ki)uj(k 2 )uo(k 3 )Y 1 5(ijo(ki) + u(k 2 ) - uo(k 3 ))5(ki + k 2 - k 3 ) 
uj{k 3 ){W{k 2 )W{k 3 ) + W{kx)W{h) ~ W{h)W{k 2 ) 

+W(k 1 )W(k 2 ) - Wik^Wih) - W(k 3 )W(k 2 )) = . (7.4) 

We used here the symmetrization of 2u(ki) to uj(ki) +u(k 2 ), the energy conservation 
us(ki) + oj{k 2 ) = u(k 3 ) in term (I), and the cyclic substitution k\ — > k 3 , k 3 — > k 2 , 
k 2 — * k\ in term (II). 

If the ergodicity condition (E) holds, energy is the only conservation law, see the 
discussion at the end of Section [T21 

(ii) Entropy. Following (|3.12|) . up to a constant, the local entropy at position r and 
time t on the kinetic scale is defined through 

s(r,t) = / dklogW(r,k,t). (7.5) 
Jt 3 

It satisfies the semi-conservation law 

^s(r,t) + V-j s (r,t)=a(r,t) (7.6) 

with the entropy flow 

j s (r, t) = (27T)- 1 / dkVuj(k) log W(r, k, t) (7.7) 
Jt 3 

and the entropy production 



a(r, i) = 7 / dkidk 2 dk 3 (u(ki)uj(k 2 )uj(k 3 )) 1 5(u(k 1 ) + u(k 2 ) - u(k 3 )) 
J -> 

5{h + k 2 - k 3 )W{r, h,t)W(r, k 2 , t)W{r, k 3 , t) 

(W(r, k^t)- 1 + W(r, k 2 , t)- 1 - W(r, k 3 , t)' 1 ) 2 . (7.8) 

Clearly a > 0. To derive the expression ()7.8|) one uses the same identities as for the 



26 



energy, 



7 / dk 1 dk 2 dk 3 (u(k 1 )uj(k 2 )uj(k 3 ))- 1 W(k 1 y 1 {25(uj(k 1 ) + u(k 2 ) - u(k 3 )) 

S(h + k 2 - k 3 ) (W{k 2 )W{k 3 ) + W{k x )W{k 3 ) - W(h)W(k 2 )) 
+5(u(ki) - uj(k 2 ) - cu(k 3 ))S(k 1 — k 2 — k 3 ) 
(W(k 2 )W(k 3 ) - WikJWfa) - WikjWfa))} 



= 7 / dkidk 2 dk 3 (uj(ki)uj{k 2 )uj(k 3 )) 1 5(cu(k 1 ) + u(k 2 ) - u(k 3 ))5(ki + k 2 - k 3 ) 
Jt 9 

W(k 1 )W(k 2 )W(k 3 )(W(k 1 )- 2 + W(k 2 )~ 2 + 2W(k 1 )~ 1 W(k 2 )- 1 

^wih^wih)- 1 + w{k 3 y 2 - w{h)- 1 w{k 2 )- 1 - wihy'wih)- 1 ) 

= a. (7.9) 

The entropy production vanishes if and only if 

Wih)- 1 + Wih)- 1 - W(h + h)- 1 = (7.10) 

on the set {(k\, k 2 ) G M. 6 \u(ki) + w(k 2 ) = u(ki + k 2 )}. As will be discussed in 
Section H21 if the ergodicity condition (E) holds, the only solution to (j7.10j) is 

(3 > is a free parameter. Physically, {3 is the inverse temperature, j3 = (k^T)^ 1 . 
We will use temperature units such that fee = 1. 

(in) Stationary solutions. For the spatially homogeneous Boltzmann equation under 
the ergodicity condition (E) the only stationary solutions are of the form (J7.ll)) . If 
there would be another stationary solution, its entropy production has to vanish, 
which means ()7.10j) has to be satisfied, in contradiction to ()7.11|) being the only 
solution of flTTUJ) . 

(17. 1 1ft is in accordance with equilibrium statistical mechanics. In thermal equi- 
librium, the nonlinearity can be ignored in the kinetic limit and the Gibbs dis- 
tribution is Z~ l exp[— (3Hq\. This is a Gaussian measure with Wigner function 

W p {k) = {(3u{k))-\ 

To have the one-parameter family (|7.11|) as the only stationary solutions is a 
remarkable prediction of the phonon Boltzmann equation. It means that the weak 
nonlinearity thermalizes the gas of phonons. For example, one could set up an initial 
state with nonvanishing phonon current, j n (0) = (27r) _1 J J3 dkVu>(k)W(k,t = 0) ^ 
0. Through umklapp processes this current degrades in the course of time and 
lim^oo j n (t) = 0. If there are no umklapp processes, as for the continuum wave 
equation below, on the kinetic level there are stationary states which maintain a 
constant phonon current. 
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8 Wave turbulence 



Wave turbulence has become a generic term for, possibly multicomponent, wave 
equations with weak nonlinearity. Examples are listed in and include waves 
on liquid surfaces, acoustic turbulence, and the nonlinear Schrodinger equation for 
dispersive media. The link to our discussion comes from the fact that apparently 
kinetic theory is the most powerful method available to handle the nonlinearities, 
a concrete field of application being the dynamics of ocean waves [TI]. The under- 
lying physical space is K 3 , possibly 1R 2 , which means that we briefly return to the 
continuum setting from the end of Section |21 For wave turbulence, typically one 
is interested in a stationary nonequilibrium state which is sustained by pumping in 
energy at large scales and dissipating it at small scales. Thus the focus is on sta- 
tionary solutions of the spatially homogeneous equation with the appropriate source 
terms added. Here we only discuss the derivation of the kinetic equation from the 
Klein-Gordon equation ([2.20|h 

1)2.20)1 has the dispersion relation u(k) = (u 2 + k 2 ) 1 ^ 2 , luq > 0. For three- wave 
interactions the resonance condition reads 

uih) + uj{k 2 ) = u{h + k 2 ) , (8.1) 

where momentum conservation, ki + k 2 = k 3 , has been used already. If u> > 0, then 

v(ki + k 2 ) < uj(ki) + u(k 2 ) (8.2) 

and (j8.1j) cannot be satisfied. If Uq — 0, the vectors must be collinear which again 
yields a vanishing collision term. Thus on the kinetic time scale we have to turn to 
four- wave interactions in order to have the nonlinearity still in effect. 
The "simplest" example is 



& 
Of 



2 



Ve\<p 3 . (8.3) 



As in ()8.2|) one concludes that the merging of three phonons into one and the splitting 
of one phonon into three are forbidden processes on the kinetic scale. The only 
remaining possibility are pair collisions, see Figure El For the formal derivation of 
the kinetic equation one proceeds as in Section El with the result 

^-W(k) + Vu(k) ■ V r W(k) = (8.4) 
at 

—\ 2 (2n)- 3 / d^k^k^k^u^ujik^ik^ujih)) 
5(uj(k) + uj(ki) - uj(k 2 ) - Lu(ks))5(k + ki — k 2 — k 3 ) 

[Wik^Wik^Wih) + W{k)(W(k 2 )W(k 3 ) - WikjWfa) - W (ki)W (k 2 ))\ . 

Here we use the standard convention for Fourier transformation in M 3 and j d 3 k is 
understood as the integration over all of 1R 3 . 
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Figure 3: A four phonon collision with number conservation. 



In a recent series of studies Nazarenko and coworkers reconsider the derivation 
of (J8.4)) from a different perspective. We explain the method in Appendix 118.31 

The kinetic equation (|8.4|) preserves number, momentum, and energy of phonons. 
This is also reflected by the formally stationary solutions 

W a p 7 (k) = (f3u + a ■ k + 7) _1 (8.5) 

with (3 > \a\ and 

7 +(/? 2 -a 2 r 1/2 (/? 2 + a 2 )>0, (8.6) 

so to have W a /3 7 (k) > 0. The solutions ()8.5j) have infinite energy because of the 
divergence at large k. Such states have not been included in our set-up. In particular, 
starting from finite energy initial data, the system cannot properly reach thermal 
equilibrium. 

The additional conservation laws are also reflected in the size dependence of the 
thermal conductance. At the high temperature side of the sample on the average 
more phonons are created than at the low temperature side. The collisions conserve 
momentum. Thus there is a laminar flow of phonons which transports energy inde- 
pendently of the size of the sample. In distinction, a real fluid has diffusive energy 
transport, since no particles are created, resp. destroyed, at the boundary. 

To turn to the issue of wave turbulence, one considers a spatially homogeneous 
situation and augments the kinetic equation ()8.4|) phenomenologically with a driving 
term as 

§-W(k,t) = C(W(t))(k) + T(k)W(k,t) , (8.7) 

where as a shorthand the collision operator is denoted by C(W). One is interested 
in the steady state W s , for which dW s /dt = 0. From the if-theorem we know that 
/ d 3 kW s (k)- 1 C(W s )(k) > 0. Therefore 

J d 3 kT{k) < 0. (8.8) 

In addition, to have energy and phonon number conservation in the steady state it 
must hold that 

J d 3 kuj(k)W s (k)T(k) = 0, J d 3 kW s (k)T(k) = 0. (8.9) 
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We imagine to have a narrow band source of energy at small k and a sink at large 
k, compare with (|8.7|) . In the intermediate regime one has to solve then 

C(W s ) = 0. (8.10) 

To be specific let us consider the wave equation with uj{k) = \k\, i.e. u = 0. 
Since u is homogeneous, and so are the collision rates, it is natural to look for a self- 
similar solution of ()8.10|) of the form W s (k) = \k\~ a . Indeed, besides the equilibrium 
values a — 0, 1, one obtains the solutions 

W s {e \k) = \k\~ 5/3 , W s (n) {k) = \k\~ 4/3 , dimension d = 3 . (8.11) 

As their equilibrium counterpart, the solutions (|8.11|) have infinite energy because 
of ultraviolet divergence. The true steady state for ()8.7|) has the power law of (j8.11|) 
only in some intermediate regime and the at large \k\ negative T supposedly ensures 
that / d 3 kW s (k) < oo. 

The physical meaning of the steady states in (J8.11|) can be understood through 
studying the flux in energy space [TU] . It turns out that W s supports a constant 
energy flux directed from small u to large u, while Wg supports a constant phonon 
number flux directed from large uo to small uj. 



9 Quantizing phonons, locally quasifree states 

The basic Hamiltonian ()2.8|) is readily quantized by regarding q x as multiplication 
operator and substituting —id/dq x for p x as acting on the Hilbert space L 2 (R, dq x ) 
attached to the site ieZ 3 . To derive the phonon Boltzmann equation it is conve- 
nient to switch immediately to the notation of second quantization and to work in 
Fourier space rather than with the spatial lattice. The one-particle Hilbert space is 
then 

t) = L 2 (T 3 ,dk) (9.1) 
out of which we construct the bosonic Fock space through 

oo 

^=0(ff n )symm. (9.2) 

n=0 

Here (f) (X, "')symm is the n-fold tensor product restricted to wave functions symmet- 
ric under permutation of labels. On JF we define a scalar Bose field with cre- 
ation/annihilation operators a(k),a(k)*, which satisfy the canonical commutation 
relations 

[a(k),a(k')} = = [a(k)* , a(k')*} , [a(k), a(k')*} = 5{k - k') . (9.3) 

Properly speaking, one has to smear a{k) to a(f) = J T3 dkf{k)a{k) with / e fj to 
have a well-defined operator on Fock space. 
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In terms of the Bose field a(k) the quantization of H from ()2.8j) results in 

H = H + V + V 4 (9.4) 

with 



Hq = dkuj (k)a(k)* a(k) , 

JT3 




\{{2uj{k 3 )Y^{a{k 3 ) +a(-k j )*) , A' = A 2 /18^ 2 . (9.5) 
i=i 

We added explicitly the stabilizing quartic term. Then H > and we may take 
the Friedrich extension to make out of H a self-adjoint operator acting on Fock 
space. Physically, the Fock vacuum Q corresponds to the ground state of H for the 
infinitely extended lattice. States ip G T thus describe local excitations away from 
the ground state. In particular, far away from the origin the particles are in their 
state of lowest energy. H Q is normalized to have ground state energy zero. Thus 
states of finite mean H -energy, {tp, Hqi/j)^ < oo, are the finite energy excitations of 
the harmonic lattice. As discussed already, for the purpose of the kinetic limit this 
set-up is of sufficient generality. 

The role of the Gaussian measures in the classical model is taken over by the 
quasifree states. They can be defined through their moments 

(a{k)*a{k')) Q = R(k,k'), 

m n 

(n^)*]!^)^ = -Werm^^, ^)}i<^<" (9-6) 

3=1 0=1 

with perm denoting the permanent of a matrix. Clearly, the positivity of the state 
(-)Q is ensured only if R > as a quadratic form. 

If R is a projection, then the state (-) Q is pure (i.e. given by a vector in J 7 ) and 
is a coherent state in the usual terminology. Let iV denote the number of phonons, 

N= I dka(k)*a(k). (9.7) 

Then 

tiR= [ dkR{k,k) = (N) Q . (9.8) 

Thus R to be of trace class is a sufficient condition for the state (-)Q to live on Fock 
space. 
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In kinetic theory the building blocks are states which are locally translation in- 
variant, stationary under the dynamics generated by H Q , and have a strictly positive 
entropy per unit volume. The obvious candidates are quasifree states characterized 
by the covariance 

(a(k)* a(k')) Q = W(k)S(k — k') , W{k) > . (9.9) 

Such a state has infinite energy and is thus outside of Fock space. The required 
mathematical framework is well studied ^7j, but will not be needed here. Instead 
we will consider a scale of states in Fock space labelled by e such that locally a state 
of the form ()9.9|) is approximated in the limit e — > 0. 

We still need to compute the entropy per unit volume of a quasifree state, com- 
pare with (|3.1ip . ()3.12|) . We choose the periodic box [1,£] 3 and a quasifree state of 
the form (|9.6jh ()9.9|) with discrete k, k 6 • • • The corresponding density 

matrix is denoted by pi and has the entropy 

St = -trpilogpe, (9.10) 

trace over Fock space, pe is of the form Z~ x exp[— X(k)a*(k)a(k)} with A = 
log ((1 + W)/W). Therefore 

S e = ((1 + W(k)) log(l + W(k)) - W(k) log W(k)) (9.11) 

which becomes 

\im r 3 S e = / dk((l + W(k))]og(l + W(k)) - W(k)]ogW(k)) . (9.12) 

We now follow the classical intuition and give ourselves a phonon distribution 
function W(r,k). Then s > 0} is a family of quasifree states with the 

property that, defining 

W £ (y, k) = {e/2f [ drie i27ry - v (a(k - er]/2ya(k + er]/2) Q ' £ (9.13) 
for y G (eZ) 3 , one has 

hmW £ ([r\ £ ,k) = W(r,k) . (9.14) 

Note that 

V / dkW £ (y,k) = (N}^ } (9.15) 



y£{s1/2) 



which means that 



(jV) Q ' £ = e~ 3 (9.16) 

under the condition ([9.14)1 . We impose (■) CI ' £ as the scale of initial states. Adopting 
the central assumption of kinetic theory, the state (-)t/ e at the long time e~ 1 t is well 
approximated by a locally quasifree state, which, as to be argued in more detail in 
the following section, results in the phonon Boltzmann equation for the quantized 
lattice vibrations. 
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10 Derivation of the phonon Boltzmann equation 
(quantized model) 



We follow the scheme of Section Inland use atomic units. The evolution equations are 
still given by (|6.3)L now interpreted as Heisenberg equations for the quantized field. 
As major difference to Section the order of the field operators must be respected. 
Let (-)t be the state at time t under the dynamics e~ tHt with initial quasifree state 
as in (|9.fjj) . ()9.13|) . The two-point function still satisfies ()fi.4j) . However, in the 
expression (j6.1H) for G(q,p,t,s) we used the commutativity of the fields to lump 
terms together, which has to be undone in the quantum context. Also, the Gaussian 
factorization (j6.12J) is to be replaced by the expectation over the locally quasifree 
state {-) s , which amounts to, for example, 

(a(ki)a(k 2 )*a(k 3 )*a(h)) s = (a(k 1 )a(k 2 )*) s (a(k 3 )*a(k i )) s 

+ (a(k 1 )a(k 3 )*) s (a(k 2 ya(k 4 )} s + (afaMh)) a (afo)Mh)*) , . (10.1) 

Note that the last term on the right vanishes by assumption. Transferred to the 
Wigner function the ordering results in 

s 3 (a(k - sr l /2Ya(k + e V /2)) s/e = W £ ( V , k, s) , 

e 3 (a(k + er}/2)a{k - er]/2)*) s/£ = 8(r}) + W £ (r], k, s) . (10.2) 

Otherwise the computation of Section |U1 can be repeated verbatim. Perhaps 
somewhat unexpected at first glance, the collision term is modified only through a 
linear term and becomes 

7 / dk 1 dk 2 {u{k)u{k 1 )uj{k 2 ))- l (25{uj{k) +u(k t ) - u(k 2 )) 

5{k + h- k 2 ) (W{r, h)W{r, k 2 ) + W{r, k)W{r, k 2 ) - W{r, k)W{r, h)) 
+5(uj(k) - u(h) - uj(k 2 ))5(k -h- k 2 ) (W(r, h)W(r, k 2 ) - W(r, k)W(r, h) 

-W(r,k)W(r,k 2 ))^ , (10.3) 
where 

W(r, k) = 1 + W(r, k) . (10.4) 

Properties of the Boltzmann equation are more readily seen by writing the colli- 
sion operator with an apparent cubic nonlinearity. This results in the conventional 



33 



form of the phonon Boltzmann equation, 

JV(r, k, t) + ^Vu(k) ■ V r W(r, k, t) (10.5) 

= 7 I dk 1 dk 2 (u(k)u(k 1 )uj(k 2 ))- 1 \2S(u(k) + uj(k 1 )-iu(k2))5(k + k 1 -k2) 

(t^(r, fc, £)#(r, fc x , t) W(r, fc 2 , t) - W{r, k, t)W{r, h, t)W(r, k 2 , t)) (I) 
+5(uj(k) - 00(h) - Lo(k 2 ))5(k -ki- k 2 ) 

(W(r, k, t)W(r, h, t)W(r, k 2} t) - W(r } k, t)W(r, k h t)W(r, k 2 , t)) } . (II) 

The Boltzmann equation ()10.5|) is one of our central results. It reduces to the 
classical phonon equation ()4.7|) through omitting the tilde in (|10.3J) . 



11 Feynman diagrams 

The iteration leading to Equation ()6.6|) suggests to develop more systematically the 
time-dependent perturbation theory. In this section we will do the first step in a 
program which needs to be completed. For simplicity let us assume an initial state 
which is translation invariant and quasifree with covariance 

(a(Jfe)) Q = 0, (a(k)a(k')) Q = , 

(a{k)*a{k')) Q = 5{k - k')W(k) , (11.1) 

compare with (|9.6|) . By the magic of Wigner functions, a slowly varying initial 
measure would require small modifications only. Since there is no spatial variation, 
kinetic scaling amounts to merely consider the long times e~ l t. By translation 
invariance 

(a(k)*a(k')) t/£ = 6(k - k')W e (k, t) . (11.2) 
As discussed already, one expects that 

limiy e (M) = W(k,t) (11.3) 

and W(k,t) to satisfy the spatially homogeneous version of the phonon Boltzmann 
equation (|10.5|) . Let us set 

W(k,l) = l + W(k), W(k,-l) = W(k) (11.4) 

and, as before, 

<j>{k, h, k 2 ) = X(8u(k)u(ki)u(k2))~ 1/2 . (11.5) 
Then the phonon Boltzmann equation (|10.5|) is written more concisely as 

d f 

—W(k,a,t) = 4ir > / dk 1 dk2(f>(k,k 1 ,k 2 ) 2 

8{au + a x uJi + a 2 uo 2 )5(ak + a t ki + a 2 k 2 ) 

(W(h, a h t)W(k 2 , a 2 , t) + aaiW(k, a, t) W ( k *> ? ' *)) ( 1L6 ) 

a=±l 
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with initial conditions from (jll.ljl . Here we use as shorthand uj = cu(k), uj(ki) = uji, 
Lj(k 2 ) = uj 2 . Note that the term with a,ai,cr 2 = 1 vanishes, since u(k) + oj(k\) + 
uj(k 2 ) > 0. 

To derive ()11.6|) from the microscopic dynamics, we take the time-dependent 
perturbation theory as starting point. The Heisenberg equations for the quantum 
field are given by ()6.3|) with the shorthand (j6.1|) . Inserting them in time-integrated 
form yields the identity 

m mm 

(JJa(fc i ,a j ))t = exp [it ( a j uj{k j ))] ( JJ a(%, a j )) Q 
j=i j=i j=i 

+iyfe / ds exp [i(t — s) ( Gju{kj))\ 

m „ 

(J2 Yl ae dk'dk"cj)(k e , k', k")5(-ak £ + o'k! + a"k") 

1=1 a ',a"=±l ^ T6 

£-1 m 

{(Y\_ a ( k j' a j)) a ( k '' a ') a ( k "' cr ") n a {kj>,vj>)s} ■ 

Note that the operator ordering is properly maintained. To generate the perturba- 
tion series for W £ (k, t) one starts with m = 2. Then on the right there is a product 
of 3 a's for which one substitutes (jll.7j) with m = 3, etc. Finally one averages 
explicitly over the initial quasifree state (-) Q . This yields 

oo 

(a(q, a q )*a(p, a p )) t/e = 5(a q , -a p )8(q - p) (W(q, a q ) + ^ W^(q, cr q , t)) . (11.8) 



n=l 



Let us postpone the issue of the convergence of the sum over n to the end of this 
section and first discuss W% for each n separately. 

S(a q /, —(J p )S(q — p)W^(q, <7i, t) is a sum of oscillating integrals. The summation 
comes from three sources 

- the sum over a', a" in (|11.7|) 

- the sum over i in (jll.7|) 

- the sum over all oriented pairings due to the average in the initial quasifree state, 



2n 



j=l pairings tt,tt' i=l 

Oriented means that in (a(fc 7r (j), o r 7r (j))a(A; 7I -'(j), cr vr /( i ))) Q on the right hand side the 
operators appear in the same order as on the left hand side. Since the integrals 
have a rather complicated structure, it is convenient to visualize them as Feynman 
diagrams. 
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Figure 4: Example of a Feynman diagram at order n = 2. 



A Feynman diagram is an oriented graph with labels. We first construct the 
graph and then the labelling. The graph uses as "backbone" In + 2 equidistant 
horizontal level lines labelled from to 2n + 1. The graph itself consists of two 
binary downward trees. The roots are two vertical bonds from line In + 1 to In. 
These bonds are continued downwards. At level m there is exactly one branch point 
with two branches. Branches do not cross, see Figure EJ Thus at level there 
are 2n + 2 vertical bonds (branches). They are connected according to the pairing 
rule resulting in n + 1 pairs. Thereby the graph consists of internal lines and two 
roots (external legs). The Feynman graph is oriented, with lines pointing either up 
(a = +1) or down (a = — 1). If there is no branching the orientation is inherited 
from the continuing vertical bond in the level above. At a pairing the orientation 
must be maintained. Thus at level a branch with an up arrow can be paired only 
with a branch with a down arrow. If the pairing is pointing to the left, it corresponds 
to the order (a*a)^ in (jll.9|) . while a pairing pointing to the right corresponds to 
(aa*)®. By construction, each internal line has two orientations and starts and ends 
at a branch point. 

Next we insert the labels. The level lines to In + 1 are labelled by times 
< t\ . . . < t2n < t. The left root carries the label q while the right root carries the 
label p. Each internal line is labelled with a wave number k. 

To each Feynman diagram one associates an integral through the following steps. 

(i) The time integration is over the simplex < t% . . . < tm < t a s dt\ . . . dt 2n . 

(ii) The wave number integration is over all internal lines as J dk\ . . . J dk K , where 
k = 3n — 1 is the number of internal lines. 

(iii) One sums over all orientations of the internal lines. 
The integrand is a product of 3 factors. 
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(iv) There is a product over all branch points. 

At each branchpoint there are a root, say wave vector k% and orientation a x , and 
two branches, say wave vectors k 2 , k 3 and orientations a 2 , cr 3 . Then each branch 
point carries the weight 

S(-a x k x + cr 2 k 2 + <T2,kz)(Ti(j){ki, k 2 , k 3 ) . (11.10) 

If one regards the wave vector k as a current with orientation a, then ()4.fi)) expresses 
Kirchhoff's rule for conservation of the current. 

(v) By construction each bond carries a time difference t m+x — t m , a wave vector k, 
and an orientation a. Then to this bond one associates the phase factor 

exp[z(t m+ i - t m )auj(k)/e\ . (H-H) 

The second factor is the product of such phase factors over all bonds. 

(vi) The third factor of the integrand is 

n+l 

llWfotTj), (11.12) 

3=1 

where k x , . . . , k n+x are the labels of the branches between level and level 1. <jj = 1 
if the pairing line is oriented to the right and <jj = —1 if oriented to the left. 

(vii) Finally there is the prefactor (— l) n e~ n . 

S(a q , —a p )S(q—p)W^(q, t) is the sum over all integrals corresponding to Feynman 
graphs with 2n+2 horizontal time slices, given the external legs q,p with orientations 
a q , a p . 

To illustrate the method let us consider the case n = 1. There are then (2 • 3) • 
3 • 4 = 72 diagrams. There is a group of 24 diagrams for which each of the two trees 
branches once. Among them there are 8 diagrams with the two trees disconnected. 
They yield the term S(q)5(p)0(e) provided u>(0) > 0. According to the rules listed, 
the remaining 16 diagrams sum up to the oscillating integral 

dkidk 2 (j)(q, fci, k 2 )(j)(p, h, k 2 ) 



It = e~ x 28{(T q , -a p ) f dt 2 I dt x V / 



CT1,(T2=±1 

(exp [i(t 2 - t 1 )(a g uj(q) - a x u{kx) - a 2 u(k 2 ))/e\ + c.c.) 

x5(-a q q + a x kx + a 2 k 2 )S(-a p p - a x kx - a 2 k 2 )W(k x , -a x )W(k 2 , -a 2 ) . (11.13) 
The limit e —>■ is covered by the argument from (|6.2(J|) and 

lim = 5(a q , — (Tp)5(q — p)Air / dt 2 V" dk x dk 2 (j)(q, k x , k 2 ) 2 

5(a q uj(q) + a x uj{k x ) + a 2 u(k 2 ))5(a q q + a x k x + a 2 k 2 ) 
W(k x ,a x )W{k 2 ,a 2 ). (11.14) 
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Secondly there is a group of 48 diagrams for which one of the two trees does not 
branch. Among them there are 16 diagrams which have an internal line with k = 0. 
They cancel amongst each other by symmetry. The remaining 32 diagrams sum up 
to I~. Its oscillating integrals are handled as for /+. Thereby one obtains 



We note the analogy with the discussion in Section |H1 The computation there is 
more lengthy, since spatial variation is included. G ga i n corresponds to I+, Gi oss to 
I~ , Gsubi to the 8 diagrams with both trees branched, and G su b2 to the 16 diagrams 
with only one tree branched. 

Kinetic theory claims that at any order diagrams divide into leading and sub- 
leading. The subleading diagrams vanish in the limit e — > while the leading ones 
have a finite limit. In fact the leading diagrams can be characterized very concisely. 

Kinetic Conjecture: In a leading Feynman diagram the Kirchhoff rule never forces 
an internal wave number i.e. a factor of the form S(kj) with some wave vector kj. 
In addition, the sum of the 2(n — m+1) phases from the bonds between level lines 2m 
and 2m + 1 vanishes for every choice of internal wave numbers. This cancellation 
must hold for m = 0, . . . , n. 

By a tricky combinatorial argument |18j it can be shown that the sum of all 
leading, according to the Kinetic Conjecture, diagrams satisfy a set of differential 
equations, which in analogy to the kinetic theory of gases is called Boltzmann hierar- 
chy. Let (/i, / 2 , . . .) be a vector of functions where f n (ki, a±, . . . , k n , a n ) is symmetric 
in its arguments. We define the collision operator C„ jn+ i through 




(11.15) 



n 




(C njn+1 f n+1 )(k 1 ,ai,...,k n ,o- n )=4ir^2 ^2 / dk'dk"4>(k £ ,k' ,k") 2 



l=\ o-',<t"=±1 



5(a e uj e + a'uj' + a"uj")5(a e k e + a'k' + a"k") 
[f n +i(h, 01, ... , ki-i, at-i, k', a',..., k", a") 




(11.16) 



5=±1 



Then the Boltzmann hierarchy reads 

~nfn(t) = Cn,n+lfn+l(t) ■ 



(11.17) 



Note that the result from (|11.14j) . (|11.15J) can be stated as 



lim(/ £ + + / £ -)=t(C 1 , 2 / 2 )(fc 1 ,a 1 ) 

£^0 



(11.18) 
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provided one sets 



f 2 (k u ax, k 2 , a 2 ) = W(h, a 1 )W(k 2 , a 2 ) . (11.19) 

The Boltzmann hierarchy has the property that initial factorization of f n is 
maintained in time, 

n 

fn{h,<7l,---,kn,Vn,t) = JJ /(%, <7j , t) (11.20) 

i=i 

and each factor satisfies the Boltzmann equation 

?-f(k,a,t) =4ttA 2 V / dk'dk"(8uu'u")- 1 

5{au + <tV + a"u")6(ak + a'k' + a"k") 
[f(k', a', t)f(k", a", t) + aa'f(k, a, t)(f(k", 1, t) + f(k", 

For the particular choice 

f(k, 1) = 1 + W(k) , f{k, -1) = W(k) (11.22) 

Equation ()11.21|) agrees with the phonon Boltzmann equation (|11.6JI . Thereby the 
Kinetic Conjecture amounts to the assertion 

limWZ(k,a,t) = ^-t n (C h2 ...C n , n+1 f n+1 )(k,a) (11.23) 

with the initial f n factorized as in ()11.19j) and single factor (jll.22)) . 

The difference between the quantized theory, discussed so far, and the classical 
theory is surprisingly minor when viewed on the level of Feynman diagrams. The 
classical a-field commutes, which however does not simplify the structure of the 
diagram. Only in the average of the initial state quasifree is replaced by Gaussian 
which according to (J11.9)) induces the modification (aa*) G = (a*a) G . Thus the 
classical phonon Boltzmann equation is obtained by setting the initial conditions 
for the hierarchy as 

n 

/„(*!, ai,...,k n ,a n )=Y[ W(kj) . (11.24) 

i=i 

One checks that (|11.21|) indeed coincides with ()4.7|) . 

So far we avoided the issue of the convergence of the series in (|11.8J) . At order 
n there are (2n — l)!(2n)!2 n /n! Feynman diagrams. If W is bounded, then a single 
Feynman diagram is of order c n t 2n /(2n)\ with some suitable constant c. Thus, 
unless cancellations are used, even at finite e the sum over n does not converge. The 
situation improves in the kinetic level. At order n there are only (48) n n! leading 



(11.21) 

-M))]- 
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diagrams. If e max of (|5.8J) is bounded, then each diagram is of order c n (e max ) n t n /n!. 
Therefore in the limit the sum over n in (|11.23|) converges provided t is sufficiently 
small. 

We conclude that the most immediate project is to establish (jll.23|) . which 
means that all subleading diagrams vanish in the limit e — > 0. This would be a 
step further when compared to the investigation [TH] , see also [23 |2T] • Of course a 
complete proof must deal with the uniform convergence of the series in (lll.8j) . 



12 Properties of the quantum phonon Boltzmann 
equation 

The Boltzmann equation (jl0.5|) for the quantized phonons differs somewhat from 
its classical cousin (|4.7|) . But the basic properties remain unaltered. As before, 
energy is locally conserved. The entropy functional has to be modified, but results 
again in a positive entropy production. Most significantly the stationary distribution 
functions are now the one-parameter family of Bose-Einstein distributions at zero 
chemical potential, 

W p (k) = (e^W - I)" 1 . (12.1) 

In the high temperature limit, (3 — > 0, they reduce to which are the 

stationary solutions of the phonon Boltzmann equation for classical lattice dynamics. 
In the sequel we discuss each item separately. 

(i) Energy. The properties (J7.1j) to (|7.3|) remain intact. One only has to show that 
j j6 dku{k)C{W){k) = 0. Inserting the collision operator from f)10.5j) one obtains 

/ dk 1 dk 2 dk 3 (uj(k 1 )uj(k 2 )uj(k 3 )y 1 {25(uj(k 1 ) + u(k 2 ) - uj(k 3 ))5(kt + k 2 - k 3 ) 

^{^(Wik^Wi^Wih) - W^W^Wiks)) + 5(u(k 1 ) - u(k 2 ) - u(k 3 )) 
8(h -k 2 - k 3 )u(k 1 )(W{k 1 )W(k 2 )W{k 3 ) - W(kt)W(k 2 )W(k 3 ))} 

dkidk 2 dk 3 (u(ki)uj(k 2 )uj(k 3 )Y l 5(uj(ki) + uj(k 2 ) - uj{k 3 ))8{ki + k 2 - k 3 ) 

u{k 3 ){W{k 1 )W{k 2 )W{k 3 ) - W(h)W(k 2 )W(k 3 ) + Wik^Wik^Wih) 
-Wik^Wik^Wih)) = , (12.2) 

where as in (17. 4j) in the first summand we symmetrized 2u(ki) to u>(ki) + u>(k 2 ) 
and used energy conservation, while in the second summand we employed the cyclic 
substitution ki — > k 3 , k 3 — > k 2 , k 2 — >• k\. 

(ii) Entropy. As can be seen from (I9.12J1 . the local entropy is defined through 

s(r,t)= dk(W(r,k,t)logW(r,k,t)-W(r,k,t)logW(r,k,t)) . (12.3) 
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It satisfies the semi-conservation law 

d 



dt s(r,t) + V ■ j s (r,t) = a(r,t) (12.4) 



with a positive entropy production a. The entropy flow is easily deduced to 

js(r 



t) = (2-kY 1 [ dkVuo{k) (W{r, k, t) log W{r, k, t) - W(r, k, t) log W{r, k, t)) . 

(12.5) 

To compute the entropy production we symmetrize as in the case of the energy, the 
role of oo(ki) being taken over by log (W(r,ki,t)/W{r,kx,t)). Then 



a(r, t) = 7 / dhdhdhicoik^cuih^ik^yHiuih) + cu(k 2 ) - co{k 3 )) (12.6) 
S(k 1 + k 2 - k 3 )f(W(r, h,t)W(r, k 2 , t)W{r, k 3 , t), W{r, h, t)W{r, k 2 , t)W{r, k 3 , t)) 
with 

f(x,y) = (x -y)\og(x/y) . (12.7) 
Clearly, <r(r, t) > unless x = y in (|12.7j) . i.e. unless 

W(kx)W(k 2 )W(k 3 ) = W(k!)W(k 2 )W(k 3 ) (12.8) 

on the set of (k\, k 2 , k 3 ) J s satisfying the constraints 

u(h) + u(k 2 ) = u>(k 3 ) , k 1 + k 2 = k 3 , (12.9) 

where we regard W and u as continued periodically to all of M 3 . Thus a = if and 
only if (|P2~%|) holds for all {k 1} k 2 , k 3 ) G R 9 on the set defined by (fT2~3jl . 
The total entropy is 

S(t) = [ d 3 rs(r,t). (12.10) 



From (ITT41 it follows that 

j t S{t)> 0, (12.11) 
which is the analogue of Boltzmann's H-theorem. 

(in) Stationary solutions. We consider a spatially homogeneous system for which 
the Boltzmann equation (jlO.Sjl reduces to 

d 

—W = C{W). (12.12) 

By definition a stationary solution has to satisfy C{W) = 0. This equality looks 
rather unapproachable and a better strategy is to use that for a solution to be 
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stationary its entropy production has to vanish. To make the resulting functional 
equations (|12.8|) and (|12.9|) more tractable we introduce as auxiliary quantity 

V> = \og(W/W) . (12.13) 

Then f|12.8|) becomes additive as 

i/>{k 1 )+if>(k 2 )=il>(k 3 ). (12.14) 

Proposition 12.1 Let the ergodicity condition (E) be satisfied and let det(Hesso;) = 
at most on a set of codimension 1. Let ip : T 3 — > R be twice continuously differen- 
tiable and satisfy the functional equation 

ipih) + il>(k 2 ) = i)(h + k 2 ) (12.15) 

on the set = {(ki, k 2 ) G T 6 | u(kx) + uj{k 2 ) = uj{k\ + k 2 )}. Then necessarily 

ip(k) = au(k) , a G R . (12.16) 

Proof: In spirit we follow Cercignani and Kremer [22J. We set k\ — k , k 2 = k'. The 
collisional invariant ip(k) -\-ip{k') is constant on the set {(k, k') G T 6 | u(k) + u(k') = 
const , k + k' = const}. Therefore there exists a function : t x T 3 -> 1 such that 

ij){k) + ip(k') = <t>(u(k) + oj(k') ,k + k'). (12.17) 
We set k = (jfe 1 , k 2 , k 3 ), u = u(k), J = cu(k'), 

d u <f>(u,k) = d<f>(u,k)/du, d a <f)(uj,k) = d<f)(uj,k)/dk a , (12.18) 

a = 1, 2, 3, and differentiate ()12.17|) with respect to k, k'. Then 

d a ip(k) = d u (p(uj + u', k + k')d a uj + d a (j)(u> + u>', k + k') , 

d a ip(k') = d u <j>(u + u',k + k')d a uj' + d a ct){uj + tu',k + k') . (12.19) 

Subtracting and symmetrizing yields 

(d a i/>(h) - d a ij){k')) {dpu(k) - d p u>(h>)) 

= (d^(k) - dpip(k')) (d a u(k) ~ d a u{k')) . (12.20) 

Differentiating with respect to k, 

d a d 1 ip(k){dj3Ui(k) — dpu(k')) + (d a ip(k) — d a ip(k')) df}d^uj(k) 

= dpd y ip(k)(d a u(k) - d a uj{k')) + (dp^{k) - d p ip(k'))d a d^uj(k) , (12.21) 

and once more differentiating with respect to k', 

d a d 1 ijj(k)d l 3dsuj(k') + d a dsil}{k')di3d 1 uo(k) 

= d p d^(k)d a d 6 u(k') + dpd s il)(k')d a d y uj(k) , (12.22) 
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which holds on A u . 

Let A = {k G T 3 , det(Hessu;(A;)) ^ 0}. As proven in Appendix [TO if k, k' G A, 
then the only solution to ()12.22j) reads 

d a d/3ip(k) = a{k)d a dpuj{k) , d a d^(k') = a(k)d a d p oj(k') (12.23) 

with some constant a(k) independent of k'. We choose now k" G A linked through 
a collision to k' and conclude that also 

d a d P i){k") = a{k)d a dpoj{k") . (12.24) 

By the ergodicity condition (E) the relation ()12.24|) extends to 

d a df3ip(k) = ad a dfsuj(k) (12.25) 

on A with some constant a. By continuity ()12.25|) extends to all of T 3 . Integrating 
(112. 25j) yields ip(k) = auo(k) + b-k + c. b = by continuity of ip and c = by (112.151) . 

□ 

Remarks: (i) Presumably our result holds under weaker smoothness assumptions on 
ip. The difficulty is that in ()12.22|) k and k! are constrained variables. 

(ii) The ergodicity condition may fail because at given k no collision is admitted 
by energy conservation. But there are more subtle cases. For example Z 3 could be 
partitioned into two sublattices which are dynamically disconnected, i.e. the elastic 
constants a(x) couple only within each sublattice. Then, at best, each sublattice 
thermalizes by itself and ergodicity is violated. 

(iii) Assume that there is some function, ip(k), such that f j3 dki[)(k)W(r, k,t) sat- 
isfies a local conservation law in the form (|7.2j) . Then the corresponding current is 
necessarily 

Mr,t) = (27V)- 1 [ dk{Vu{k))^{k)W{r,k,t) (12.26) 
and it must hold that 

/ dkip(k)C(W)(k) — (12.27) 
for all W. Repeating the computation in ()12.2j) leads to 

dk 1 dk 2 dk 3 (u 1 0J 2 ^3)~ 1 {*P(h) +^{k 2 ) - ip(k 3 ))5(uj l + u 2 — u 3 ) 
5(kt + k 2 - k 3 ) (W{k 1 )W{k 2 )W{k 3 ) - Wik^Wik^Wih)) = . (12.28) 

Hence ip is a collisional invariant in the sense of Proposition 112.11 Under the as- 
sumptions stated there, it follows that ip(k) = au(k) and energy is the only local 
conservation law. 
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For the case at hand, ui{k) > and W(k) > 0, which implies a < 0. Thus we 
have shown that under the ergodicity condition (E) the only stationary solutions of 
the spatially homogeneous Boltzmann equation are 

W p {k) = (e*"W - (3>0. (12.29) 

f3 is fixed through the initial condition as L 3 dku(k)W(k, t = 0) = L 3 dku(k)Wp(k) 
by conservation of energy. Thermodynamically (3 is the inverse temperature and the 
entropy of Wp according to ()12.3|) is the equilibrium entropy of an ideal Bose gas at 
zero chemical potential. 



13 The linearized collision operator 

The thermal conductivity, in the kinetic limit, is determined through the inverse 
of the linearized collision operator. We will explain the standard argument in the 
following section. Here we merely study the linearized collision operator as a linear 
operator in L 2 (T 3 ,dk). 

We consider the spatially homogeneous Boltzmann equation ()10.5|) . which we 
write as 

jW = C(W). (13.1) 

Under the ergodicity condition (E) the only stationary solutions of (j!3-H) are the 
thermal Wp. We fix (3 and linearize at Wp, where the convenient way of writing the 
perturbation is 

W = Wp + WpWpf . (13.2) 
To linear order in /, (113. lj) then becomes 

WpWp^f = -Lf (13.3) 
with the linearized collision operator 

Lf(k) = -7 f dk l dk 2 {u{k)u{k l )uj{k2)Y 1 (13.4) 

(28(u(k) + u(h) - u(k 2 ))5(k + h- k 2 ) {(Wp(k 2 ) - Wp(h))Wp(k)Wp(k)f(k) 

+(Wp(k 2 ) - Wpik^Wpik^Wpik^fih) + (Wp(k) + Wp(h))Wp(k 2 )Wp(k 2 )f(k 2 )) 
+5(u(k) - u(ki) - u(k 2 ))5(k -h-k 2 )(- (Wp(h) + Wp(k 2 ))Wp(k)Wp(k)f(k) 

+{Wp{k 2 ) - Wp{k))Wp{k^W p {k x )f{k 1 ) + (Wp(k 1 ) - Wp(k))Wp(k 2 )Wp(k 2 )f(k 2 ))y 
In each term we use the 5-constraint which leads to identities of the type 

Wp(h)Wp(k 2 )Wp(h) = Wp{kx)Wp{k 2 )Wp{h) on u{ki) + oo(k 2 ) = tu(k 3 ) . (13.5) 
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Then (I13.4j) simplifies to 



Lf(k) = 7 / dfadk 2 (u(k)u(fa)u(k 2 ))- L 25(u(k) + u{fa) - tu(k 2 ))5(k + fa - k 2 ) 

W f3 (k)Wp(fa)W f} (k 2 ) (f{k) + f{fa) - f(k 2 )) + 5{u{k) - u(fa) - uj(k 2 )) 
5(k -fa- k 2 )Wp(k)W p (fa)Wp(k 2 ) (f(k) - f(fa) - f(k 2 ))) . (13.6) 



Let (■,■) denote the inner product in L 2 (T 3 ,dk). Using once more (J13.5)) . the 
quadratic form for L is given by 

= 7 / dfadk 2 dfa(u(fa)u)(k 2 )u)(fa)y l 

5{io{fa) + u(k 2 ) - u(k 3 ))S(fa + k 2 - fa)W p (fa)Wp(k 2 )Wp(fa) 
(g(fa) + g(k 2 ) - g(k 3 )) (f(fa) + f(k 2 ) - f(k 3 )) . (13.7) 

Thus it is evident that 

L*=L, L>0, Lto = 0. (13.8) 

Note that any zero eigenvector of L, (f, Lf) = 0, must be a collisional invariant 
in the sense of ()12.15j) . Hence, under the stated assumptions, the eigenvalue zero 
is nondegenerate. In the classical limit, (5 — > 0, Wp and Wp are to be replaced by 
(/5c<j) _1 . Then L equals the linearization of 1)4. 7|) . 

The spectral properties of L have not been studied, to our knowledge. But they 
seem to fall into the standard folklore picture of kinetic theory. L can be written as 



Lf(k) = - / dk'A(k, k')f(k') + V(k)f(k) . (13.9) 
The "potential" follows from (J13.6)) as 

V(k) = ^Wpitycuik)- 1 [ dfa{oj{fa)u{k + fa))' 1 (2S(u{k) + co(fa) - tu(k + fa)) 
Wp(fa)Wp(k + fa) + 5(co(k) - uo{fa) - uo(k + fa))Wp(fa)Wp(k + fci)) , (13.10) 
while the integral kernel A has the form 

A(k, k') = (13.11) 
2 7 { - {uj(k)uj(k')uj(k + k'))- l Wp{k)Wp{k')Wp{k + k')8(u(k) + u(k') - u(k + k')) 
+ {u(k)uj(k')uj(k - k'^W^Wp^W^k - k')5{uj{k) - u(k') + u(k - k')) 
+ (u(k)u(k')u(k - k'^W^kW^k'W^k - k')5(u(k) - u(k') - u{k - fa))} . 

Under our assumptions on u the potential is bounded away from zero, < c_ < V(k) 
but not bounded sine uj = 0. For k fixed, A(k, k') is concentrated on a set of 
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codimension 1. The kernel of A 2 is a function, but A 2 (k, k') has singular points, in 
particular A 2 {k, k) = oo. We conjecture that trA 4 < oo. If so, the bottom of the 
continuous spectrum of L is c_. L has a spectral gap and the continuous spectrum 
extends to oo. On the linearized level the homogeneous system relaxes exponentially 
fast to equilibrium. 



14 Thermal conductivity 

We look for a stationary solution of the Boltzmann equation 1)10.5)1 . to say 

{2ir)- 1 Vu-V r W = C(W), (14.1) 

which has approximately a linear temperature profile T(r) = + VT • r with 
|VT| <C 1. Of course, the ergodicity condition (E) has to be imposed. On the left 
hand side in ()14.1j) we assume local equilibrium in the form (e w ^ k ^ T ^> — 1) _1 while 
on the right hand side we expand W = Wp + WpWpf . Then ()14.1j) becomes 

(Vu ■ VT)uW p W^ 2 = -Lf . (14.2) 

Since, as argued before, the zero eigenvalue of L is nondegenerate and the corre- 
sponding eigenvector u is orthogonal to (^u)ujWpWp oo, L can be inverted and 

/ = -(27i)- 1 f3 2 L- 1 W f3 W l3 uj(VT ■ Vcj) . (14.3) 
The steady state heat (=energy) flux is then 

j e . = (27T)- 1 / dk{W p + WpW p f)uVu 

= (2ti)- 2 P 2 (ujVujW (3 W( 3 , L^WpWpUjVuJ ■ VT) . (14.4) 
The thermal conductivity k is defined through Fourier's law j e = — kVT, hence 

Kaa'(T) = p 2 {2Ti)- 2 {W p Wpu\7 a uj } L- l WpWpu\7 al u) , /3 = 1/T. (14.5) 
For the case at hand, k is diagonal, K aa i = S aa 'K, and 

«(T) = \T- 2 (2ix)- 2 (WpWpuVuo , L- x W p WpU ■ Vuo) . (14.6) 

Inserting Fourier's law into the local conservation of energy (|7.2j) yields a non- 
linear diffusion equation for the energy transport, 

|e(r,t) = V- («(T(e))^Ve(r,t)) (14.7) 

with the thermodynamic relation 

e(T) = / dkuWp, (3 = l/T. (14.8) 
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It would be of interest to establish ()14.7jl as the hydro dynamic limit of the Boltzmann 
equation (|1(J.5|) . 

We discuss the qualitative temperature dependence of the thermal conductivity 
At high temperatures, Wp and Wp are replaced by (Pu)^ 1 . Then 

= \ T ~\2ixy 2 (-VuJ, (Ld)- 1 - ■ Voo) (14.9) 
with the classical linearized collision operator 

L c \f{k) = 7 / dkidk 2 dk 3 (u(k)ui(ki)uj(k 2 ))~ 2 

(26{u(k) + u}(ki) - u(k 2 ))8{k + h- k 2 ) (f{k) + f(h) - f(k 2 )) 
+5(co(k) - - u(k 2 ))8{k -h- k 2 ) (f(k) - f(h) - /(fc 2 ))) .(14.10) 

Therefore the temperature dependence is multiplicative and 

<T) = |, T large, (14.11) 

with 6*h determined by ()14.9|1 . 

At low temperatures the temperature dependence of k is not so easily accessible. 
For P — > oo 

W p (k) = e-^ k \ (14.12) 

which means that the number of energy carrying phonons is greatly reduced. On 
the other hand, normal processes conserve momentum and thus do not degrade the 
phonon current. Only in umklapp processes, momentum is transferred to the lattice. 
But umklapp becomes rare at low temperatures. It is argued in Chapter 2.2, 
that the latter effect dominates resulting in the exponential increase 

k(T) ~ e 0l/T , 0j > , T small . (14.13) 

It would be of interest to have bounds based directly on (J14.5)) which confirm such 
a low temperature behavior. 

For real materials the dependence ()14.13|) is not so easily resolved, since the con- 
ductivity is dominated by scattering from isotope mass disorder, as will be discussed 
in the next section. For mass purified samples the conductance is limited by the size 
of the probe. 

15 Isotope disorder 

At low temperatures the thermal conductivity is limited by impurities. Even for 
a chemically pure crystal, in their natural abundance the crystal atoms come as a 
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random mixture of isotopes. Artifically enriched, resp. purified, samples have also 
been manufactured so to provide a test of the predictions by the theory. In the kinetic 
limit, the effects of impurities and small anharmonicities are additive. Therefore 
we study here random isotope substitution in the harmonic approximation. If m x 
denotes the mass of the atom at site x, in the frame of our toy model the equations 
of motion read 

d , . 1 

-tlIxW = — Px(t) 
at m x 

-^Px(t) = - ^2 a(y - x)q y (t) - u%q x (t) , x G Z 3 , (15.1) 

compare with (J2.11j) . For isotope disorder the mass ratio is of order 10~ 2 . Therefore, 
in a good approximation we may set 

— = (l + v^) 2 , £<1, (15.2) 
m x 

with {£ x , x G Z 3 } a collection of independent, identically distributed random vari- 
ables. Let us denote by E the expectation with respect to the £ x 's, i.e. the disorder 
average. We assume E^) = and \£ x \ < c so that m x > for sufficiently small e 
as required for mechanical stability. 

To derive the kinetic equation we first follow the scheme devised for the weak 
nonlinear ity, also to emphasize that the structure is in parallel. To mathematically 
justify the decoupling step a distinct strategy is required, however, see Section [TBI 

We rewrite the equations of motion ()15.1)) in terms of the a-field as defined in 
(|2.16p . which means in terms of the homogeneous system. Since the evolution is 
linear, there is no difference between the classical and quantum model, possibly 
except for the choice of the initial state and thus the initial Wigner function. One 
obtains 



d 
It 



a(k,a,t) = iauj(k)a(k,t) — i\feo \ . \ dk\dk2{uj{k)uj{ki)) 



1/2 



x5(-ak + a-th - k 2 )(7 1 a{k 1 , a u t)£(k 2 ) + 0{e) , (15.3) 



where we take the quantum framework, to be definite. Compared to (J6.3|) in essence 
one of the a-factors has been replaced be £. Since a(k, a, t) depends on the disorder, 
the equations of motion are, so to speak, nonlinear in the couple (a, £). 
The object of interest is the Wigner function 

W £ ( V , k, t) = e 3 E((a(k - 67 1 /2)*a(k + e V /2)) t/E ) (15.4) 

on the kinetic time scale e~ l t. In (115. 4|) there are two averages, one over the disorder, 
E, and one over the initial state. To be physically consistent we think of a scale of ini- 
tial states as explained in Sections El and El Since the equations of motion are linear, 
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there is however a much wider choice. In the classical model the initial configuration 
could be deterministic. Quantum mechanically the initial wave function could be in 
the one-particle space f). All what is required is that the Wigner function (|15.4|) (pos- 
sibly substituting e 3 by some other prefactor) has a limit at the initial time t = 0. 
The disorder average is taken only to avoid extra difficulties in the derivation. Phys- 
ically one expects (a(k — erj /2)* a(k + er] / 2)) t / e to be self-averaging in the limit e — > 0. 
More precisely, the random variable J dr)dkf(r), k){a(k — er]/2)*a(k + erj/2)) t / e , f a 
smooth e-independent test function, will tend with probability one to a deterministic 
limit as e — > 0. No disorder average should be needed, in fact. 

Let us see how the arguments from Sections El and El transcribe to the present 
situation. As before the atomic scale is used. Then 

^K((a(p)*a(q)) t ) = i(u(p) - Lu(q))E((a(p)*a(q)) t ) 



+e I dsG(q,p,t — s, s) (15.5) 
Jo 



with 



G(q,p, t, s) = — E 2_, I dk\dk2 2_, / dlidl2{oj{ki)u{a\ki — k 2 )) 



1/2 



{u(h)u{nh - k)) 1/2 ak 2 )ai 2 ) (e u ^ + ^ k ^5{-p + a,h - k 2 ) 

x {8(q + nk - h^ini^h, o-i)a(/i, n)) s + S^axh + nh - 1 2 ) 

xn (a(k, n)a(q)) a ) + e u ^ + ^ k ^5(q + o x k x - k 2 ) (5(-p + n h - h) 
xcriri(a(Zi,ri)a(A;i,cri)) a + 8{-<j\k\ + nh - Z 2 )ri(a(p)*o(/i,ri)) s )J (15.6) 

The homogeneous term vanishes, since E(£ x ) = 0. 

There is no need to carry the argument any further, since we have seen it already. 
The analogue of the assumption of local stationarity is to factorize, on the kinetic 
scale, the disorder average as 

E(f(A;i)f(A:2)(a(Zi)*a(g)) 8 ) = E(e(An)f(A; 2 ))E«a(Zi)*a(g)) a ) , (15.7) 

for example. The rapidly oscillating time integral generates the 5-function for the 
energy conservation and makes terms as (aa) and (a*a*) to vanish. After these steps 
only four terms are left which combine into the linear Boltzmann equation for the 
limit Wigner function W, 

^W(r, k, t) + ^uj(k) ■ V r W(r, k, t) 

= 2nE{il)oj{k) 2 [ dk!8{uj{k)-u{ki))(W{r,ki,t)-W{r,k,t)) . (15.8) 



T 3 
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If one wants to compute the thermal conductivity including the isotope disorder, 
one merely has to add to L in (|14.2|) the impurity scattering in the form 

L-J{k) = 27rE(£ ) WWpWp [ dk x 5{u{k) - u(kx))(f(k) - f(ki)) . (15.9) 

By energy conservation L; randomizes on each energy shell. Thus the zero eigen- 
vectors of Li are of the form h(u(k)) with arbitrary h. The Planck distribution is 
singled out by the anharmonicities and can be thought of as a specific initial con- 
dition in the current context. Following the arguments in Sectional the thermal 
conductivity is given through 

k(T) = -(3 2 {2tt)- 2 (W i3 W p ujVuj, (L + L^WpWpu ■ Vu) . (15.10) 
3 

In the limit of vanishing anharmonicity, 7 — > 0, k can be computed more explic- 
itly, since, by symmetry k ~» —k, Vcu is an eigenfunction of L\. Then 

«,(T) = b^Y'C-ml))- 1 I dkW,W,{Vujf^— (15.11) 
6 2 J T 3 t{uj) 

with 

t{u)= / dh8(u - w(Jfei)) . (15.12) 

If ujq > 0, K{ vanishes exponentially as e _/3a;o . If u = 0, the competition between 
the divergence of t{uj) for small ui and the suppression of phonons results in a de- 
pendence as «i(T) — T^ 1 for T — ► 0. 



16 Mapping to a Schrodinger-like equation with 
a weak random potential 

The linear evolution equation (|15.1|) suggests to use time-dependent perturbation 
theory. Let us set 

!). "-(!*)■ (161) 



Then 

*Vp) = (A + ^'(p) (16 ' 2) 



d 

and 



00 „ 

e (A +v ^y) 4 = + ^ e"/ 2 / dt n . . . dt^-^Ve^"-'"-^ . . . Ve Ah . 

n=1 J o<ti<...<t„<t 

(16.3) 
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We insert the propagator (jl6.3j) into the definition of the Wigner function and 
average over disorder. The leading term is exp[— t(l/27r)Vc<j • V r ] on the kinetic 
scale. The term of order y/e vanishes because E(^) = and the term of order e 
yields, when kinetically scaled and taking the limit e — > 0, 

f d ae -(*-»)(l/*r)Va».V rjLie -«(l/^)Va,.V rW ( lg 4) 

Jo 

with W the initial Wigner function and Li the linear collision operator of f)15.8j) . 
Thus we only have to study systematically the higher orders of the perturbation 
series and to convince ourselves that they yield the corresponding time-dependent 
perturbation series for ()15.8j) . Unfortunately, while the principle is correct, it will 
never lead to a proof, since there are too many terms in the perturbation series. 
Even if we postulate that the {^j are independent Gaussians, the number of pair- 
ings is n\/2 n l 2 {n/2)\ which are balanced by a factor a n t n l 2 / \n/2)\ from the time 
integrations. Thus the series converges only for \t\ < t with a suitable to on the 
kinetic time scale. 

Erdos and Yau [21] study the one-particle Schrodinger equation with a random 
potential which has a mathematical structure comparable to (jl6.2j) . Thus the prob- 
lem of an exploding number of terms in the perturbation series also arises. To 
circumvent this blockage, they expand only up to iV = N(e) and estimate the re- 
mainder by using the unitarity of the unexpanded Schrodinger evolution. To copy 
their method we have to exploit that the energy 

H = \ E (c 1 + ^fpi + E °(f - x )w*) ( 16 - 5 ) 

is conserved for each realization of the disorder. The energy depends on £. This 
is rather inconvenient and we transform to new fields such that the flat £ 2 - n orm is 
conserved. Let us regard 

n xy = [ dke i2nHx ~ v) uj{k) (16.6) 

as a linear operator in £2 = ii^T?)- Under our assumption Q has an exponential 
decay in \x — y\, which possibly worsens as u — > 0. We define 



with components 



v = (v> + ,T/n (16.7) 

^ = -j={nq x ±i(i + • (16.8) 

Note that ||V' + || 2 = IIV' - !! 2 = H. Thus the £ 2 -norm of ip is conserved in time. The 
"0-field evolves according 

.±(^\ = (n o\(i>+\ r(tn + nz ^-^U^ (169) 
dtU-) vo -n U-J +ve l-en + ^ -^n-ndU-h 
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where £ x is regarded as a multiplication operator, (£f) x = £ x f x - We use the short 
hand 

i^ = (H + y/eV)i>, H e = H + ^V (16.10) 

and regard ()16.9|) as an evolution equation in £ 2 © (-I- Clearly, H £ is bounded and 
H £ = H*. Thus e~ lHet is unitary. Physical initial data are constrained to satisfy 
= ip~ , but this will be imposed only at the very end. 
Since ip is a 2-spinor, the Wigner function becomes a 2 x 2 matrix. Inserting the 
kinetic scaling, one has 

W* a ,{y,k,t) = 2- 3 [ dr]e i27Ty - r, ^ a (k-eri/2,t/ey^ a '(k + er]/2,t/e), (16.11) 

a = ±, a' = ±, with k G T 3 and y G (eZ) 3 . Note that, because of the definition 
()16.8|) . we deviated slightly from previous conventions. In particular f j3 dkW^_ + 
(y, k, t) now acquires the meaning of an energy density at kinetic time t. The off- 
diagonal element, W+_(t), picks up the fastly oscillating phase factor exp[±2iu(k)t / e] 
Hence it vanishes upon time averaging. For example, W+_(t) determines the dif- 
ference between kinetic and potential energy, which is indeed a fast variable. By 
symmetry Wt_{t) is obtained from W ++ (t) by substituting uj by — uj. Thus we only 
have to deal with W+ + (t). 

The Wigner function at fixed t typically oscillates on small scales and only upon 
integrating against a smooth test function one expects to have a limit. Thus let 
J : IR 3 x T 3 — > R be a smooth, rapidly decreasing function with its Fourier transform 
with respect to the spatial argument denoted by J. The Wigner function integrated 
against J becomes then 



(<W;+M> = / drj / dk^ + (k-E7]/2)*J(7],k)ij + (k + er]/2). (16.12) 

We now choose a sequence of initial conditions ip £ such that \\ip e \\ < const and such 
that the initial Wigner function has a limit, 

\im(J,Wl + {ip £ ]) = f J(r,k)no(drdk) . (16.13) 

In addition one has to impose tightness in the sense that 

lim limsup V |^| 2 = 0. (16.14) 

|a;|>R/e 

At this level of generality fio(drdk) is a positive, bounded measure. With ip e as 
initial datum the time-evolved field is given by 

ijj(t) = e~ iHet ip £ . (16.15) 
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Clearly, the issue is to determine (J, W E [ip{t/e)]) in the limit e — > 0. 

To achieve the existence of the limit one needs two conditions on the dispersion 
relation u. 

(i) The first condition we have met already in Section 5 and requires u to be a Morse 
function, meaning that all critical points of to are isolated and nondegenerate (no 
zero eigenvalue in the quadratic approximation). 

(ii) The second condition is the crossing estimate, which refers to the decay estimate 
of a particular oscillatory integral over T 3 x T 3 . It is too technical to be stated 
explicitly here. The crossing estimate is verified for a few particular dispersion 
relations [2H1 E3 I2H] • It is not excluded that with improved technology the crossing 
estimate can be reduced to the Morse property. 

Before stating our result we have to explain what we mean by solution of the 
Boltzmann equation (|15.8|) with a measure as initial condition. The standard 
method is to switch to the dual equation and to prove that it is a contraction 
semigroup on C(R 3 x T 3 , R), the space of bounded and continuous functions, which 
follows from the key observation that, since u is Morse, the operator 

Bf(k) = 2rrE(^)uj(k) 2 [ 5(u(k) - u{k'))f{k')dk' (16.16) 
satisfies \\Bf\\ < c\\f\\ in C(T 3 ) for some c > 0. In particular the total collision rate 

v(k) = 2ncu(k) 2 [ 5(uj(k) - uo(k'))dk' (16.17) 

is continuous, thus bounded. To the Boltzmann equation there is associated the 
stochastic process (r(t), k(t)) with state space R 3 x T 3 . It is governed by 

j f r(t) = ^Vw(fc(t)) , (16.18) 

where k(t) is a Markov jump process on T 3 with jump rate 2irK(^Q)u!(k) 2 5(u(k) — 
u){k'))dk' . We define the measure fi t (drdk) as the joint distribution of (r(t),k(t)) 
when started with /x as initial measure. 

The following theorem is the main result of a joint paper with J. Lukkarinen 

ESI 

Theorem 16.1 Let u be a Morse function and satisfy the crossing estimate and let 
ip £ G £2 © h be uniformly bounded and such that \16.1ty . \16.1$ hold. Then 

\imE((J,W £ (ij(t/e)))) = [ fx t (drdk)J(r,k), (16.19) 
where fi t is the solution of the Boltzmann equation hl5.£fy with initial datum (Xq. 
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Remark: The dispersion relations (|5.3j) and ()5.4jl with u>o > are Morse and satisfy 
the crossing estimate. 

Any description of the methods used in the proof would lead us too far astray 
Let me only emphasize that they are based on techniques developed by Erdos and 
Yau (21], see also for estimating Feynman diagrams and for cutting delicately 
the perturbation series in an e-dependent way. In a recent paper Chen considers 
the Schrodinger equation on a lattice with nearest neighbor hopping and a random 
potential V with V(x), x e Z 3 , a collection of independent random variables. His 
estimates greatly helped in our proof. Chen j57j also shows that for his model the 
convergence of the Wigner function holds in probability which is a strong indication 
that the same property should hold for isotope disrder. 

17 Guide to the literature 
17.1 Phonon Boltzmann equation 

I am not an expert in phonon physics and the guide reflects merely my own read- 
ing. The focus is deliberately somewhat narrow and I deal only with the rigorous 
derivation and a few basic properties of the phonon Boltzmann equation. 

The seminal paper on the subject is R. Peierls P from 1929. He is the first 
one to write down the phonon Boltzmann equation ()10.5|) . Nordheim [28J follows 
a similar path for weakly interacting quantum gases. Peierls' derivation consists 
in a careful application of Fermi's golden rule. His argument, with variations and 
modernized notation, has been repeated many times. A standard reference is the 
Handbuch article by Leibfried j2J. An excellent textbook discussion is Callaway [5]. 
I very much enjoyed the monograph by V.L. Gurevich He also applies the Fermi 
golden rule but in addition discusses extensively the physical conditions required for 
its applicability in the derivation of the Boltzmann equation. As a standard, the 
Fermi golden rule is introduced in the context of the spatially homogeneous system. 
Spatial variation is simply added in the most obvious way. A great advantage of the 
Wigner function formulation is to incorporate spatial variation from the outset. 

Since the most interesting aspects of phonon physics are related to quantiza- 
tion, the classical anharmonic crystal tends to be ignored. But in his basic article 
Peierls also treats the classical system. Brout and Prigogine j2H] provide a more de- 
tailed account, which is summarized in the book by I. Prigogine on nonequilibrium 
statistical mechanics [20]. He and Peierls argue that, through a random phase ap- 
proximation, the joint distribution of the a(k)*a(k) satisfies a diffusion equation in 
the high-dimensional phase space. Reducing to the one-particle distribution yields 
a nonlinear evolution equation for W(k,t), in spirit similar to the structure one has 
in the Kac model of kinetic theory E21- At the time such reasoning was very 
fashionable. But its underlying assumptions are rather dubious, see Appendix 118.31 

As regards to derivation from the microscopic Hamiltonian model the next level 
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is to improve on the Fermi golden rule, which started with the work of van Hove 
and lead into the development of diagrammatic expansions in parallel with 
similar techniques in quantum field theory. This is a vast area, still active today. A 
very readable account with focus on weak coupling and Boltzmann type transport 
equations is the slim monograph by S. Fujita [H|. He discusses the impurity problem 
and electron-electron collisions. But he could have treated phonons, as well. Fujita 
immediatly employs Wigner functions as a matter of fact, which makes one wonder 
who originally pushed this concept as a tool for transport equations. In his famous 
paper Wigner [33] introduces the notion but then applies it to the semiclassical limit 
of the quantum statistical partition function. 

The importance of local stationarity has been stressed mostly in the quarters of 
mathematical physics, since it is one central property which needs to be established 
when proving the validity of a macroscopic equation. Erdos, Salmhofer, and Yau [16J 
discuss the strongly related problem of electron collisions in the same spirit as done 
here. Benedetto et al. [TH] are more ambitious and, ignoring the issue of absolute 
convergence, study the dominant terms of the time- dependent perturbation series 
in the kinetic limit. Most likely, their techniques extend to the present case. 

The harmonic lattice with isotope disorder is in its structure rather similar to 
a one-particle Schrodinger equation with a weak random potential. We refer to 
j2H 1213 EI3 123 ESI for recent advances and the derivation of the corresponding 
kinetic equation. Our Theorem 14.1 relies on their work. 

Bal, Komorowski, and Ryzhik study the continuum wave equation with a 
weakly disordered index of refraction. They consider a high frequency approximation 
and prove that in this limit the Wigner function is governed by (|15.8|) with the jump 
collision operator replaced by a spherical Laplacian, which turns out to by the small 
angle approximation to the collision operator in ()15.8|) . 

Compared to its famous sister the phonon Boltzmann equation has received little 
mathematical attention, for no good reason. While we expect that much of the tech- 
nology developed in the context of the Boltzmann equation carries over, we point 
out that the phonon case has two simplifying features: The wave vector space is 
compact and more importantly the velocity, Vo;(fc), is uniformly bounded. How far 
this will carry, only a detailed study can show. The derivation of hydrodynamics 
should be more accessible, since the Boltzmann equation has only a single conser- 
vation law and its corresponding nonlinear diffusion equation (|14.7|) has a global 
solution, say in a finite macroscopic box with initial data bounded away from 0. 

17.2 Energy transport in anharmonic chains 

Classical anharmonic chains are a challenging test ground for the numerical inte- 
gration of Newton's equation of motion ever since the seminal work of Fermi, Pasta 
and Ulam j37j. With increasing computer power steady state current transport for 
chain lengths up to 10 4 , in exceptional cases even 10 5 , are reported. These studies 
mostly investigate strong anharmonicities and are thus only loosely related to the 
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kinetic theory discussed here. For this reason we merely refer to a few review arti- 
cles. Jackson j3H] summarizes to work up to 1978, an authorative 2003 update being 
Lepri, Livi, and Politi [39 . Bonetto, Lebowitz and Rey-Bellet j3U] emphasize more 
theoretical aspects, in particular large deviations and the fluctuation theorem. 

The numerical simulations available provide no clear evidence, whether kinetic 
theory is applicable in dimension one (and two). In the kinetic theory of gases colli- 
sions in one dimension are degenerate, since particles just pass through each other. 
On the other hand for lattice dynamics three phonon processes are non-degenerate, 
as can be checked explicitly for the dispersion relation (|5.4|) . Therefore, in general, 
the collision operator does not vanish. Ergodicity is more questionable. For our 
standard example ()5.4|) at u> = the components [— ~,0] and [0, |] are not linked 
through collisions. As Uq increases these components shrink. In their steady state 
the phonon current would not vanish. Other couplings, four phonon processes, or 
thermal boundary drive could restore ergodicity. Whether the microscopic model for 
small but fixed anharmonicity has regular energy transport remains to be studied. 
Only some loosely related results are available. Aoki and Kusnezov [IT] numerically 
simulate the case ujq = 0, A = 1 and report good evidence for normal heat con- 
duction, i.e. a steady state energy current proportional to 1/N with N the chain 
length. For the same model the current-current momentum and energy correlation 
functions are studied in [32] • A variety of other harmonic nearest neighbor chains 
with anharmonic on-site potential is investigated in |4*5] . Lefevre and Schenkel [4*4"] 
attempt to expand directly the steady state probability distribution under thermal 
boundary conditions. They report the term of order A. From the perspective of 
kinetic theory the term of order A 2 would be related to the chain length of order 
A- 2 . 

18 Appendix 

The proofs given below are due to J. Lukkarinen. 

18.1 Three phonon processes in case of nearest neighbor 
coupling only 

For nearest neighbor coupling the dispersion relation reads 

3 

u(k)= (cu 2 + 2^(l -cos(2vrF))) 1/2 . (18.1) 
i=i 

We prove that 

u(k) + uj(q) -u(k + q) > ujq/2 (18.2) 

for all q, k G T 3 . Therefore in this case three phonon collisions are prohibited by 
energy conservation. 
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We set z = (z 1 , z 2 , z 3 ) and 



z j (k) = - -L e - t27TkJ ) , a > 1 . (18.3) 



Then by direct computation \z(k)\ = u(k) with u determined uniquely by a. We 
find 

\u{k + q)- u{k)\ < \\z{k + q)\ - \z{k)\\ < \z(k + q) - z(k)\ (18.4) 



and 



a 

with z = z at a = 1. Therefore 



\z(k + q)- z{k)\ 2 = -\z (q)\ 2 < \z (q)\ 2 (18.5) 



uj(k) +uj(q) - u(k + q) > uj(q) - \uj(k + q) — oj{q)\ 

> u(q) - \z {q)\ = (co 2 + \z {q)\ 2 ) l/2 - \z Q {q)\ > co /2 (18.6) 



for all \z (q)\ 2 > 0. 



18.2 Entropy as the logarithm of phase space volume 



In (22] Garrido, Goldstein, and Lebowitz argue that whenever a suitable set of 
"macrovariables" evolves in time according to an autonomous deterministic law, 
then the entropy functional, defined as the logarithm of the phase volume associated 
to specified values of the macrovariables, is increasing in time. A system of weakly 
interacting phonons should be no exception and we will explain why. 

Notationally it is convenient to choose the wave number torus as T 3 = [0, l] 3 . 
If the lattice volume is [1, . . . ,£] 3 , then the wave numbers are discretized as k e 
(T^) 3 = (^ _1 [1, • • • ,^]) 3 - We partition the unit torus into cubes Aj of side length 5, 
5M = 1, j = 1, . . . , M 3 . Accordingly we set 

Hj = a(k)*a(k) (18.7) 

as a function on phase space (M 6 )^ 3 . The H/s are the macrovariables. They are 
assumed to take a value close to £ 3 ej with 

e- = 5 3 [ d 3 kW{k) . (18.8) 

Let e = (ei, . . . , ejv/s). The corresponding region in phase space is 
T(e, 6, v) = {(q,p) e (M 6 )'Vfe ~v)< H 3 < t\e 3 + v) , j = 1, . . . , M 3 } . (18.9) 
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Then, using the equivalence between mirocanonical and canonical ensemble, 

A-/ 3 

lim lim r 3 log|r(e,5,z/)| = (<5 3 Vloge,-) + logvr + 1 . (18.10) 

u^O I-+00 — ' 

3=1 

If one now refines the partitioning into cubes by letting 5 —>■ 0, one arrives at the 
entropy functional 

/ d 3 k(logW{k) + log7r + 1) (18.11) 
Jt 3 

in accordance with (|3.12|) . 

The quantum case is rather similar, once it is realized that the operators from 
(|18.7|) are a family of commuting operators. The conditions in ()18.9j) define a pro- 
jection operator P(e, 5, v) on bosonic Fock space and 

M 3 

lim lim £~ 3 logtrP(e, 5, v) = {5 3 ((1 + e 5 ) log(l + e ) - e 3 - logej . (18.12) 

As before, upon refining the partition by letting 8 — > one arrives at the entropy 
functional 

/ d 3 k((l + W(k))log(l + W(k)) - W{k)\ogW{k)) (18.13) 

JT 3 

in accordance with 1)9.11)1 . 

As to be expected on general grounds 03] and as confirmed by Propositon 112.11 
the thus defined entropy is increasing in time when W(k, t) evolves according to the 
phonon Boltzmann equation. 

18.3 The Brout-Prigogine equation 

In the context of wave turbulence, over recent years the validity of the Boltzmann 
transport equation has been scrutinized with the aim to understand the necessity 
for corrections [J7J HH1 EH!- One part of the enterprise are numerical simulations 
testing the validity of Gaussian local statistics [2011211. In these works the authors 
follow the Brout-Prigogine scheme [2E1 ED] • Since it differs from our approach, to 
comment on their method might be instructive. 

We consider the finite volume A = [1,...,£] 3 C Z 3 . With periodic boundary 
conditions our Hamiltonian reads 

H =\Hti + \ E " P (^)« + ^5>L (18.14) 

where a p are the periodized elastic constants and includes Uq. We now rotate 
q, and p, such that a p (x,y) becomes diagonal. It has the eigenvalues oj\ with 
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k G A* = {t 1 [1, . . . , £]) 3 , the dual lattice. If q k , Pk denotes the new coordinates and 
momenta, we further switch canonically to action-angle variables through 

qk = (Jk/uJ k ) 1/2 cos a k , Pk = {u) k J k ) 1/2 sin a k , (18.15) 
< Jki ctfc G 27rT. In action-angle variables the Hamiltonian becomes 

H = ^Jk + y/iH^a) . (18.16) 

The precise form of H\ is easily worked out, but not needed for our summary. The 
equations of motion are then 

d 

a k = u k + y/e—H^J, a) , 

OJk 

j k = -V^^-H 1 (J,a). (18.17) 

da k 

Clearly, the u;'s are the fast variables while the actions change slowly. 

We impose the initial distribution, p(0), on phase space which evolves under 
the flow f)18.17|) to p(t). p(0) is taken to depend only on J, the random phase 
approximation, and one is interested in the distribution of slow variables at the 
kinetic time e~ x t, 



/•2tt 

(J) = HiW 1 / dxkjpUa^-H). (18.18) 



Brout and Prigogine use second order perturbation theory for the Liouvillean, which 
suggests that pf (J) evolves approximately by a diffusion process. The computation 
is very readably explained in [28, pp. 36-60] and there is no need to repeat. As net 
result they obtain a diffusion process on (1R+)^ 3 with generator 

Lf{J)=l 22 ( UJ kUk'U k »y 1 d~(uJk + uj k > - LO k f>)5(k + k' - k") 

k,k',k"£A* 

r JL + JL _ JL\ Jk j k ,j k jJL + JL _ JL)fU) . (i8.i9) 

K dJ k dJy dJk>> Wfc dJ k > dJk»' 

(118. 19j) is a system of coupled diffusion processes. If we consider one representa- 
tive triple, (Ji, J2, J3) G (R+) 3 , then the diffusion process (Ji(t), Jiit), J%{t)) moves 
along the line {(Ji(0), J 2 (0), J 3 (0))+A(1, 1, -1), A G R}. The diffusion process never 
exits the domain (M+) 3 , since the diffusion coefficient, J1J2J3 vanishes sufficiently 
fast towards the boundary. 

According to (|18.19|) the first moment evolves as 

^-(Jk)t = (LJ k ) t . (18.20) 
dt 
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Taking I — > oo and assuming the factorization (Jk'Jk")t — {Jk')t(Jk")t one arrives 
at a closed equation for (Jk)t- As a check on consistency, it indeed agrees with the 
Boltzmann transport equation (J4.7|) . 

The tricky part of the argument is the diffusion approximation (jl8.19J) . For fixed 
£, the limit e —>■ in ()18.17|) is covered by the perturbation theory for integrable 
systems, see e.g. Chapter 5. Even if the initial phases are assumed to be random, 
there is simply no diffusion approximation in sight. The motion of the angles is 
quasi-periodic, thus much too regular for the purpose of diffusion. One is forced 
to take with e — > simultaneously £ — > oo. Kinetic scaling requires £ = O^ 1 ), 
which means to enter murky waters. It remains to be seen whether there is some 
intermediate scale on which ()18.19|) is a valid approximation. 

To my understanding, the transformation to action-angle variables easily misses 
the central physical mechanism for the validity of the kinetic description. It is 
the wave propagation in physical space, and its good spatial mixing properties, 
which ensures that even in presence of a small nonlinearity the wave field retains 
approximately the Gaussian statistics. 

18.4 Solutions to ( 112.221 ) 

We set 

d a d/3ip(k) = A a p , d a d/3ip(k') = A a p , d a dpu(k) = B a p , d a dpu(k') = B a p . (18.21) 

Then (I12.22j) transcribes to 

A ai Bp & + A aS Bp y = ApyBas + Ap S B aj (18.22) 

and we have to find out all possible solutions under the condition that B and B are 
invertible. We multiply in (|18.22|) with (i? _1 ) 7 y and (S -1 )^/ and sum over 7, 5. 
Let us define 

C = AB-\ C = AB- 1 . (18.23) 
Changing 7', 8' back to 7, S yields 

Ca-ySps + CasSpy = Cj3yS a s + Cps5 ay . (18.24) 

In f|18.24|) we choose indices a ^ (3 ^ 7, where it is used that d > 3, and we 
set 5 = a, resp. 5 = (3. Then C a p = c a 5 a p. Correspondingly from a 7^ (3 7^ 5 and 
7 = a, resp. 7 = j3, it follows that C a p = c a 5 a p. Thus 

c a 5 a7 8ps + c a S a s5p 7 = cpSp^Sas + cpSpsS ai . (18.25) 

Setting a — (3 — 7 = S one concludes c a = c a and setting a = (3, = 7, (3 = 5 one 
concludes c a = cp, a 7^ (3. Combining both identities, there exist some constant a 
such that 

C a p = a5 a p , C a p = a5 a p (18.26) 
and consequently, using (|18.23|) . 

A = aB, A = aB. (18.27) 
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